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ABSTRACT 

We continue our series of papers on open cluster distances with a critical assessment of the accuracy of main- 
sequence fitting using isochrones that employ empirical corrections to the color-temperature relations. We use 
four nearby open clusters with multicolor photometry and accurate metallicities and present a new metallicity 
for Praesepe ([Fe/H] = +0.1 1 ±0.03) from high-resolution spectra. The internal precision of distance estimates 
is about a factor of 5 better than the case without the color calibrations. After taking into account all major 
systematic errors, we obtain distances accurate to about 2% - 3% when there exists a good metallicity esti- 
mate. Metallicities accurate to better than 0. 1 dex may be obtained from BVIcK s photometry alone. We also 
derive a helium abundance for the Pleiades of Y = 0.279 ± 0.015, which is equal within the errors to the Sun's 
initial helium abundance and that of the Hyades. Our best estimates of distances are (m-M)o = 6.33 ±0.04, 
8.03 ±0.04, and 9.61 ±0.03 to Praesepe, NGC 2516, and M67, respectively. Our Pleiades distance at the 
spectroscopic metallicity, (m— M)q = 5.66 ± 0.01 (internal) ±0.05 (systematic), is in excellent agreement with 
several geometric distance measurements. We have made calibrated isochrones for -0.3 < [Fe/H] < +0.2 
available online. 

Subject headings: Hertzsprung-Russell diagram — open clusters and associations: individual (M67, 
NGC 2516, Pleiades, Praesepe) — stars: distances — stars: abundances — stars: evolu- 
tion — stars: activity 



1. INTRODUCTION 

The determination of accurate distances is the key to under- 
standing how stars and the Galaxy have formed and evolved. 
From protostars in star-forming regions to ancient tracers of 
the halo, improved distances have refined stellar evolution- 
ary theory and Galactic structure models (e.g., Reid 1999). 
The Hipparcos mission (Perryman et al. 1997) was especially 
valuable, providing trigonometric parallaxes for ~ 10 s stars 
to precision of 1-2 mas (Perryman et al. 1997). These par- 
allaxes, however, are only useful for individual stars within 
^100 pc. Most open clusters are much more distant than this 
"horizon," but a half-dozen of the nearest clusters have 10-50 
or more Hipparcos stars, yielding cluster parallaxes ostensi- 
bly accurate to 5% or better (Mermilliod et al. 1997; Perryman 
et al. 1998; Robichon et al. 1999; van Leeuwen 1999). 

Main-sequence (MS) fitting, also known as the photomet- 
ric parallax method (e.g., Johnson 1957; Siegel et al. 2002), 
has long been used to estimate distances to individual stars 
and star clusters beyond the limits of parallax studies, and is 
considered to be a robust and well-understood technique. It 
was therefore a big surprise when the Hipparcos distances to 
the Pleiades and Coma Ber open clusters were in disagree- 
ment with distances from the MS-fitting at more than a 3 a 
level (Pinsonneault et al. 1998). It is difficult to reconcile a 
short Pleiades distance with stellar interior and spectroscopic 
abundance studies. A high helium abundance would make 
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a cluster fainter than expected from its metallicity, and solu- 
tions of this type have been discussed in the literature for the 
Pleiades (Belikov et al. 1998). However, this is difficult to 
understand since there do not seem to be nearby field stars of 
similar characteristics in the Hipparcos catalog (Soderblom 
et al. 1998), and the helium enhancement would have to be 
enormous (Y 0.34). In addition, it has been suggested that 
the metal abundance from spectroscopy may have been sig- 
nificantly overestimated (Percival et al. 2003). An argument 
was also made that distance estimates from theoretical stel- 
lar models have been overestimated for young clusters due to 
unknown, age-related physics (van Leeuwen 1999). 

However, the most likely explanation is related to the Hip- 
parcos parallaxes themselves. Pinsonneault et al. (1998) 
showed that the 12 bright stars near the center of the Pleiades 
all had virtually the same parallax, ~ 9 mas, more than 1 mas 
larger than the mean parallax for other cluster stars. They at- 
tributed this to a local zero-point error of the individual stellar 
parallaxes that are correlated over the Hipparcos' 0.°9 field 
of view (van Leeuwen & Evans 1998). These quasi-random 
errors were caused by the Hipparcos great-circle data reduc- 
tions, as Makarov (2002, 2003) proved by re-reducing the 
Pleiades and Comar Ber cluster parallaxes in a different way 
that correctly obtains the absolute zero point of parallax. Ad- 
ditional effects may result from the way the Hipparcos data 
were obtained and analyzed, and a more elaborate reduction 
of the Hipparcos parallaxes promises to produce improved 
distances and better understood errors (van Leeuwen 2005; 
van Leeuwen & Fantino 2005). 

The discrepant Hipparcos result for the Pleiades subse- 
quently led to many efforts to determine the cluster's distance 
from binaries and independent parallax measurements (e.g., 
Munari et al. 2004; Pan et al. 2004; Johns-Krull & Ander- 
son 2005; Soderblom et al. 2005). These results support the 
longer distance scale from MS fitting, verifying that the Hip- 
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parcos result was in error. With a formal error of ~ 1 % from 
these measurements, the Pleiades represents a second system 
(besides the Hyades) with a sufficiently accurate distance for 
a precision test of stellar evolutionary models. 

Even though the controversy over the Pleiades distance is 
now settled, a critical assessment of the MS fitting technique 
is still required to reliably estimate a distance. In fact, MS fit- 
ting using theoretical isochrones is a complex process that in- 
volves both physical and empirical considerations (e.g., Stauf- 
fer 2001). There are, however, many opportunities to check 
the construction of the isochrones. Stellar evolution models 
can be tested against the Sun and other stars, such as eclipsing 
binaries, that have accurate masses and radii. Furthermore, 
multicolor photometry in nearby clusters and field stars can 
be used to test the bolometric corrections and color-effective 
temperature (7^) relations to transform theoretical quantities 
(luminosity and T^f) to magnitudes and colors (e.g., Vanden- 
Berg & Clem 2003). 

In our first two papers of this series (Pinsonneault et al. 
2003, 2004, hereafter Paper I and Paper II, respectively), we 
began a long-term effort to assess the accuracy of distances 
from MS fitting and to reduce or eliminate systematic er- 
rors in the process, particularly those arising from the trans- 
formation of theoretical to observational quantities. In Pa- 
per I, we demonstrated that stellar models from the Yale Ro- 
tating Evolutionary Code (YREC; Sills et al. 2000) are in 
good agreement with masses and luminosities for the well- 
studied Hyades eclipsing binary vB 22 (Torres & Ribas 2002). 
These models also satisfy stringent tests from helioseismol- 
ogy, and predict solar neutrino fluxes in line with observa- 
tions (Basu et al. 2000; Bahcall et al. 2001; Bahcall & Pin- 
sonneault 2004). In Paper II, we showed that the models pro- 
vide a good match to the spectroscopically determined tem- 
peratures (Paulson et al. 2003) for individual Hyades mem- 
bers with good parallaxes (de Bruijne et al. 2001). However, 
we found that any of the widely-used color- r ef f relations (e.g., 
Alonso et al. 1995, 1996; Lejeune et al. 1997, 1998) fail to 
reproduce the observed shapes of the MS in the Hyades; dif- 
ferences in broadband colors were as large as ~ 0.1 mag. The 
existence of these systematic errors in the colors in the pres- 
ence of agreement between the spectroscopic and theoretical 
L— r e ff scales strongly implies that there are problems with the 
adopted color- T e n relations instead of errors in the theoretical 
r e ff scale. Therefore, we proposed empirical corrections to the 
color- r e ff relations from Lejeune et al. (1997, 1998) that were 
adopted in the isochrone computations. 

In this study, we generate a set of isochrones over a wide 
range of age and metallicity, and test the validity of the 
Hyades-based color-r e ff corrections using extensive multi- 
color photometry of four well-studied nearby open clusters. 
We show that isochrones employing the Hyades empirical 
corrections precisely match the observed MS shapes, except 
where anomalously blue colors in young open clusters have 
been previously noted (Stauffer et al. 2003). Furthermore, we 
demonstrate that the empirical corrections improve the inter- 
nal precision of the isochrones by examining the consistency 
of distances derived from several color-magnitude diagrams 
(CMDs). 

We also assess various sources of systematic errors in the 
MS-fitting technique. Previously, Pinsonneault et al. (1998) 
considered the effects of age, metal abundance, helium, red- 
dening, and systematic errors in the photometry, demonstrat- 
ing that these could not explain the short distance to the 
Pleiades from Hipparcos. Terndrup et al. (2002) paid atten- 



tion to the adopted reddening law in a discussion of the dis- 
tance to NGC 2516. Here we extend the error analysis more 
quantitatively, emphasizing photometric calibration issues 
and the bias in distance estimates induced by the presence 
of unresolved cluster binaries or field foreground/background 
stars. 

This paper also explores the effect of metallicity on the 
luminosity of the MS. Metallicity changes isochrone lumi- 
nosities more strongly than many other input parameters, and 
the degree of sensitivity depends on the color index used. 
This permits a purely photometric derivation of the metal- 
licity (e.g., Pinsonneault et al. 1998, 2000; Stello & Nissen 
2001; Terndrup et al. 2002), which can be compared to metal- 
licities derived from high-resolution spectra. An agreement 
between the photometric and spectroscopic metallicities, as 
we find in this paper, provides supporting evidence that the 
effects of metallicity on the theoretical quantities (L, T e ff) and 
on the color- T e ff relations are correctly computed. 

The distances in this paper are tied to the Hyades distance 
at (m-M) = 3.33 ±0.01 (d = 46.34 ±0.27 pc), the cluster's 
center-of-mass inferred from the Hipparcos catalog (Perry- 
man et al. 1998). Unlike the controversial Hipparcos distance 
to the Pleiades, the large angular diameter of the Hyades on 
the sky makes the cluster parallax less vulnerable to the spatial 
correlation of the Hipparcos parallax (Narayanan & Gould 
1999a,b; de Bruijne et al. 2001). 

In § 2 we compile cluster photometry, metallicities, red- 
dening estimates, and information on binarity and member- 
ship, and present a metallicity for Praesepe from new high- 
resolution spectra. In § 3 we briefly describe the construction 
of the isochrones. In § 4 we compute the distances to the 
sample clusters with the reddening fixed at previously known 
values and demonstrate that the empirical corrections improve 
the internal precision of the isochrones. In § 5 we simultane- 
ously solve for the cluster metallicity, reddening, and distance 
from the x 2 minimization. In § 6 we evaluate the effects of 
several systematic error sources, including those from cluster 
binaries and field star contamination. In § 7 we discuss sev- 
eral implications of our results. In the Appendix we address 
issues on the photometric zero points of the empirical Hyades 
isochrone. 

2. OPEN CLUSTER DATA 
2.1. Selection of Clusters 

We consider four nearby Galactic open clusters with exten- 
sive multicolor photometry: Praesepe (= M44; NGC 2632), 
the Pleiades (= M45), M67 (= NGC 2682), and NGC 2516. 
The choice of these clusters was motivated by several factors. 
All have well-determined estimates of metal abundance and 
reddening against which we can compare photometrically- 
derived values. Samples in Praesepe, the Pleiades, and M67 
are dominated by known cluster members, so we can exam- 
ine whether the Hyades-based color calibration from Paper II 
generates isochrones that precisely match the shapes of the 
MS in these clusters. Praesepe has extensive information on 
binarity, so systematic errors in distances arising from the 
presence of binaries can be explored. M67 and NGC 2516 
each have modern photometry from two independent studies, 
from which we gauge the sizes of errors that arise from pho- 
tometric calibration issues. In addition, NGC 2516 has a rel- 
atively high reddening compared to the other clusters, which 
allows us to evaluate the consequences of adopting particular 
reddening laws. 

The Pleiades is a special case because its distance has re- 
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TABLE 1 

Recent Measurements of the Pleiades Distance 



Reference 



Method 



(m-M) 



Narayanan & Gould (1999b) 
Gatewood et al. (2000) 
Makarov (2002) 
Munari et al. (2004) 
Pan et al. (2004) 
Zwahlen et al. (2004) 
Johns-Krull & Anderson (2005) 
Soderblom et al. (2005) 
Southworth et al. (2005) 
Weighted mean 



Moving cluster 5.58 ±0.18 

Ground-parallax 5.58 ±0.12 

Hipparcos reanalysis 5.55 ±0.06 

Eclipsing binary (HD 23642) 5.60 ± 0.03 

Astrometric binary (Atlas) 5.65 ± 0.03 

Astrometric binary (Atlas) 5 . 60 ± . 07 

HST parallax 5.66 ±0.06 

HST parallax 5.65 ±0.05 

Eclipsing binary (HD 23642) 5.72 ± 0.05 
5.63 ±0.02 



cently been accurately measured from astrometric and eclips- 
ing binary studies and from ground- and space-based paral- 
laxes, allowing a precise test of distances derived from MS 
fitting. Individual measurements of these studies are summa- 
rized in Table 1, and the weighted average distance from these 
measurements is (m-M)o = 5.632 ±0.017. The Hipparcos 
distances to the Pleiades and other clusters are discussed in 
§7.2. 

We use cluster ages from Meynet et al. (1993) for the 
Pleiades (100 Myr), M67 (4 Gyr), and NGC 2516 (140 Myr). 
Praesepe is generally considered to be the same age as the 
Hyades (e.g., Mermilliod 1981b); in Paper II, we assumed an 
age of 550 Myr for the Hyades, as derived from isochrones 
without overshooting (Perryman et al. 1998). As we demon- 
strate, the cluster distances are insensitive to the choice of age. 

2.2. Photometry 
2.2.1. Praesepe and the Pleiades 

We compiled optical photometry for Praesepe and the 
Pleiades mainly from WEBDA (Mermilliod & Paunzen 
2003) 5 and the Open Cluster Database. 6 Data in V and 
B-V for Praesepe came from Johnson (1952), Dickens et al. 
(1968), and Castelaz et al. (1991). Following the suggestion 
by Dickens et al. (1968), we added +0.002 to V and +0.006 
to their B-V to match the Johnson (1952) data. We adopted 
the photometry from Castelaz et al. (1991) without alteration; 
the average difference in B-V is +0.002 ± 0.012, Johnson 
(1952) being redder, among 39 stars in common (including 
nonmembers). We compared the numerous sources of V and 
B-V in the Pleiades, but did not find statistically signifi- 
cant differences in any sample. At B-V > 0.8, the scatter in 
measurements from several sources often exceeded the stated 
photometric errors; this may result from brightness and color 
changes from stellar spots on rapidly rotating stars in young 
open clusters (Stauffer et al. 2003). For stars in this color 
range, as everywhere else, we simply averaged the several 
available colors and magnitudes. 

The situation in V-Ic is less straightforward because data 
for Praesepe and the Pleiades, like those for the Hyades, came 
from either the Johnson (V - If) or the Kron (V - Ik) sys- 
tems. The colors on the Johnson system were transformed 
to the Cousins (V - Ic) system using an updated transforma- 
tion equation as described in the Appendix. The Kron colors 
were transformed to the Cousins system using the cubic poly- 
nomial derived by Bessell & Weis (1987). 

5 See http://obswww.unige.ch/webda/webda.html. 

6 See http://spider.ipac.caltech.edu/staff/stauffer/opencl/index.html. 



In Praesepe, stars with V < 12 have V -I photometry on 
the Johnson system from Mendoza (1967) and Castelaz et al. 
(1991), while fainter stars have photometry on the Kron sys- 
tem from Upgren et al. (1979), Weis (1981), Stauffer (1982b), 
and Mermilliod et al. (1990). Intercomparisons showed that 
the Mendoza (1967) and Castelaz et al. (1991) colors were on 
the same system. The same is true for the Kron data except 
that the photometry from Upgren et al. (1979) required a red- 
ward shift of +0 . 03 in V - Ik to match the other photometry. In 
the Pleiades, stars with V < 10 have V -I photometry on the 
Johnson system from Johnson et al. (1966), Mendoza (1967), 
and Iriarte (1969), while fainter stars have photometry on the 
Kron system from Stauffer (1982a), Stauffer (1982c), Stauf- 
fer (1984), Stauffer & Hartmann (1987), Stauffer et al. (1989), 
Prosser et al. (1991), and Stauffer et al. (1994). Direct com- 
parison between the Mendoza (1967) photometry and that of 
Iriarte (1969) showed that they agree well for V-Ij < 0.5, 
but the Mendoza (1967) data are systematically redder by 
0.02 mag for stars with V— Ij > 0.5. We opted to adjust the 
red Mendoza (1967) photometry to place them on the Iriarte 
(1969) scale. 

2.2.2. M67 and NGC 2516 

For M67, we used BVIc photometry from Montgomery 
et al. (1993, hereafter MMJ93) and also from Sandquist 
(2004, hereafter S04); these are analyzed separately. S04 pro- 
vided a comparison between the two, revealing statistically 
significant differences between the two studies. In Table 2 we 
compile the mean differences in the photometry for M67 and 
for other clusters as is discussed below. The first column of 
the table lists the cluster name, and the sense of the compar- 
ison is shown in the second column. The last three columns 
display the mean difference and its standard error in V, B-V, 
and V - Ic, respectively. In comparison to S04, the MMJ93 
data are fainter in V and are redder in V -Ic- As shown in 
Figure 5 of S04, the differences are largest for those stars with 
V-/ c <1.0. 

For NGC 2516, we have independent photometry in BVIc 
from Jeffries et al. (2001, hereafter JTH01) and from Sung 
et al. (2002, hereafter S02). Neither study compares their pho- 
tometry with the other. In Figure 1 we plot the differences in 
the photometry from the two studies against right ascension, 
which reveals a position-dependent difference in V (but not in 
the colors). The sense of the difference is that the photome- 
try towards the east is fainter in the S02 study. According to 
S02, photometry in this portion of the sky was obtained on a 
non-photometric night which was then adjusted to match the 
data in the rest of their survey using stars also observed under 
good conditions. As their paper lists only the average values, 
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TABLE 2 
Differences in the Photometry 



Cluster Comparison <AV) (A(B-V)) (A(V-/) C ) 



M67 MMJ93-S04 +0.017 ±0.002 +0.009 ±0.002 +0.022 ±0.001 

NGC2516 S02-JTH01 +0.016 ±0.001 -0.003 ±0.001 +0.011 ±0.001 

Hyades Ground - Tycho-l" +0.012 ±0.002 -0.009 ±0.001 

Praesepe Ground - Tycho-l a +0.009 ±0.007 +0.004 ±0.003 

Pleiades Ground - Tycho-l a -0.017 ±0.003 +0.000 ±0.002 

r.m.s. ■■■ 0.015 0.006 0.017 



"Computed for stars with Vt <9. 
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FIG. 1. — Comparison of NGC 2516 photometry as a function of right 
ascension (in degrees). The differences are in the sense of the S02 minus the 
JTH01 photometry. 

we chose to use the photometry only from regions obtained 
entirely under photometric conditions. However, there still 
remain significant differences, as shown in Figure 2, between 
the JTH01 and the S02 photometry even when the eastern data 
in the latter study are removed. Compared to the JTH01 val- 
ues, the S02 data are brighter in V and redder in V—Ic at the 
top of the MS, but are fainter in V and bluer in B-V at the 
faint end. The mean differences are summarized in Table 2. 

2.2.3. Assignment of Random Errors 

Most of the collated photometry lacks errors for individual 
stars, and those errors reported by MMJ93, JTH01, and S02 
are typically from a small number of repeat observations. Our 
MS-fitting procedure (§ 4.1) first removes stars that are sta- 
tistically separated from the MS, and this in turn requires a 
suitable error for each star. 

In Praesepe and the Pleiades, we divided the data into bins 
in V, each 2 mag wide, and then computed the median of 
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FIG. 2. — Comparison of NGC 2516 photometry, excluding stars observed 
under non-photometric conditions by S02. The differences are in the sense of 
the S02 minus the JTH01 photometry. 

rms in magnitude and colors for the relatively few stars with 
multiple measurements. We assigned this median value to all 
data points in each bin as random photometric errors. For 
NGC 2516, we followed the same binning procedure, then 
removed a systematic trend in the differences between the 
JTH01 and S02 studies by subtracting a linear function in V. 
We then computed the rms of the differences in V , B - V, and 
V -Ic in each bin, divided these by \/2, and assigned these 
errors to both data sets assuming that each study would have 
about the same errors. For M67, we followed the same pro- 
cedure, first matching the stars by their coordinates, and com- 
puting the rms of the differences after removing a systematic 
trend in V. S04 reports the error for each star based on a large 
number of measurements. We used those errors directly in 
our analysis of the S04 data, and assumed that the scatter be- 
tween the two studies was dominated by measurement errors 
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TABLE 3 

Differences in the Tycho Photometry 



Cluster 


Comparison (ABt) 


(AV T ) 


Hyades 

Praesepe 

Pleiades 


Tycho-2 - Tycho- 1 -0.009 ± 0.002 
Tycho-2 -Tycho-1 -0.018 ±0.004 
Tycho-2 -Tycho-1 -0.014 ±0.003 


-0.006 ± 0.002 
-0.025 ±0.004 
-0.016 ±0.003 


Note. — 


Computed for stars with Vt < 9. 





in MMJ93. We set the errors for the MMJ93 study from the 
rms of the differences with respect to the S04 data. 

2.2.4. 2MASSK S 

We calculated V -K s colors from the All Sky Data Release 
of the Two Micron All Sky Survey (2MASS) Point Source 
Catalog (PSC). 7 Here K s designates the "short"-/T filter in 
2MASS (Carpenter 2001). Based on PSC flag parameters, 
we excluded stars that were saturated or undetected. We also 
ignored blended or contaminated sources. The V-K s errors 
were taken as the quadrature sum of V errors and the catalog's 
"total" photometric uncertainties in K s . 

2.3. Systematic Errors in the Photometry 

Our analysis of systematic errors in the MS-fitting method 
(§ 6.1) shows that calibration errors in the photometry con- 
tribute significantly to the overall error budget. In addi- 
tion to the direct comparisons between the studies of M67 
and NGC 2516, we can use the Tycho-1 photometry (van 
Leeuwen et al. 1997) to check whether photometry of the 
Hyades, Praesepe, and the Pleiades is on a consistent system 
(the other clusters are too distant to have many stars in the 
Tycho catalog). Here we assume that the Tycho photometry 
is consistently on the same scale in all parts of the sky. In 
Table 2, we display the mean differences between the ground- 
based photometry for these clusters and the Tycho values. The 
latter were transformed to the Johnson system using the equa- 
tions in Oja & Evans (1998) for V and in Mamajek et al. 
(2002) for B-V. The values in the table were computed for 
stars with Vj < 9. 

The comparison between the Tycho-1 values and the 
ground-based data is shown in Figure 3. The scatter is usu- 
ally dominated by errors in the ground-based photometry for 
V < 8 and by the Tycho photometry for fainter stars. The 
solid line in each panel denotes equality, while the dashed 
line shows the mean difference from Table 2. 

In Table 3 we list the differences between the Tycho-1 and 
Tycho-2 data (H0g et al. 2000) for bright stars in the Hyades, 
Praesepe, and the Pleiades. The differences were computed 
directly with the Tycho values (again for W < 9), without 
transformations to the Johnson system. Taken together, the 
data in Tables 2 and 3 show that the various studies probably 
could have calibration errors on the order of 0.01 - 0.02 mag 
in V, but smaller in B-V. The situation in V-Ic is less well 
determined, mainly because we have fewer studies to com- 
pare. 

In Paper II we also noted that the Tycho-1 values for V 
were brighter than the ground-based data by +0.012; we opted 
to add this constant to the Tycho data before averaging with 
the ground-based photometry. Had we chosen to adopt the 
space-based scale instead, our calibrated isochrones would be 

7 See http://www.ipac.caltech.edu/2mass/. 



TABLE 4 

Spectroscopy in Praesepe 





Tat 












ID a 


(K) 


(cm s~ 2 ) 


(km s" 1 ) 


loge 




[Fe/H] 


KW23 


5700 


4.5 


0.5 


7.78 


0.05 


0.10 


KW58 


5850 


4.5 


0.5 


7.80 


0.06 


0.12 


KW 304 


5625 


4.4 


0.4 


7.80 


0.06 


0.12 


KW336 


5600 


4.5 


0.4 


7.79 


0.06 


0.11 



a KW = Klein Wassink (1927). 



brighter in V and also bluer in V — K s since the V and K s data 
were obtained independently. This correction was not applied 
to the isochrones used in this paper. In the Appendix we dis- 
cuss other issues involved in the isochrone calibration. 

2.4. Membership and Binarity 

We also collated information on binarity and membership in 
Praesepe and the Pleiades from WEBDA and the Open Clus- 
ter Database. Any star that was listed as a nonmember for 
any reason (e.g., from photometry or radial velocities) was 
rejected. In Praesepe, there are extensive data on binarity 
(e.g., Mermilliod & Mayor 1999; Bouvier et al. 2001). Stars 
that were designated as likely or probable binary stars were 
flagged. In the Pleiades, a considerable number of stars have 
binarity and membership information (Jones 1981; Stauffer 
et al. 1991; Mermilliod et al. 1992; Schilbach et al. 1995; 
Bouvier et al. 1997; Belikov et al. 1998; Moraux et al. 2001). 
We rejected any star with a membership probability less than 
50% but kept any star that did not have a membership prob- 
ability. S04 provides a list of single cluster members based 
on proper motion surveys and known binary systems in the 
literature and further elimination of unrecognized binaries in 
CMDs from his high-precision photometry. We used this in- 
formation for both MMJ93 and S04 photometry in this paper, 
unless otherwise specified. 

2.5. Metallicity and Reddening 
2.5.1. New Spectroscopy in Praesepe 

Obtaining accurate distances relative to our calibrating clus- 
ter, the Hyades, requires accurate relative metallicities be- 
cause the luminosity of the MS is sensitive to the metal abun- 
dance. As part of this study, we determined a new metallicity 
for Praesepe using the same method used to find the Hyades 
abundance (Paulson et al. 2003). The relative abundance of 
Praesepe with respect to the Hyades should therefore be accu- 
rate since systematic errors arising from the solar abundances, 
employed model atmospheres, oscillator strengths, or effec- 
tive temperature scales would be minimized. 

We obtained spectra of four Praesepe stars with the Magel- 
lan Inamori Kyocera Echelle (MIKE) spectrograph (Bernstein 
et al. 2003) on the Magellan 6.5 m Clay telescope at Las Cam- 
panas Observatory. A 0.35"wide slit gave a resolving power 
of <~ 55, 000 per resolution element (4 pixels) and wavelength 
coverage from 4500 to 9200 A. The spectra are of high qual- 
ity, with S/N > 100. 

The spectra were reduced using standard IRAF 8 packages. 
The stellar parameters - T e ff, surface gravity (logg), microtur- 

8 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 



6 



An et al. 



O 
CO 

a. 

m 
I 

x) 
G 
3 
o 

> 



> 

I 

PQ 
< 



0.2 



0.0 



-0.2 



0.2 



0.0 



-0.2 



1 1 1 1 1 1 1 1 

o 

- Hyades o - 

° - 


^ 1 , , 1 , , 1 , ^ 

- Praesepe °o- 


1 1 1 1 1 1 1 c 
- Pleiades °- 


° o 

" 1 , , 1 , , 1 , - 


° oo _ 



o 

1 , , 1 , , 1 , 


o 00 o ° 

° p °t 

S Q 



— O— 
1 1 1 1 1 1 1 1 o 



_ 1 1 1 1 1 1 1 1 _ 

- Hyades 

CP 


_ 1 1 1 1 1 1 1 1 _ 

- Praesepe 

o _ 

: ' 


1 1 1 1 1 1 1 1 

— _Q 
© 

- Pleiades s 

o- 

°o °" 


o-o^e 8^889#^^« - - 

O 

1 1 1 1 1 1 1 1 ° 


o °o- 


° -C 

o 
o 

" 1 , , 1 , , 1 , ~ 


" 1 , , 1 , , 1 , " 



TABLE 5 
Cluster Metallicity from 
High-Resolution Spectroscopy 



V 

ground 

FIG. 3. — Differences in V (top) and B — V (bottom) between the ground-based and the Tycho-1 photometry, after transforming the latter into the Johnson 
system (see text). The differences are in the sense of the ground minus the Tycho-1 values. The solid line indicates equality, while the dashed line is plotted at 
the weighted mean difference for stars with Vr < 9 (Table 2). The scatter for V > 8 mainly reflects errors in the Tycho photometry. 

bulence (£), and [Fe/H] - were derived following the proce- 
dure in Paulson et al. (2003), and are listed in Table 4. We 
employed the spectral synthesis code MOOG (Sneden 1973) 
and used stellar model atmospheres based on the 1995 ver- 
sion of the ATLAS9 code (Castelli et al. 1997). Within IRAF, 
we used Gaussian fits to 15 Fe I lines and nine Fe II lines, a 
subsample of those listed in Table 1 of Paulson et al. (2003). 
The T e ff was derived by requiring that individual line abun- 
dances be independent of excitation potential and that micro- 
turbulence (£) be independent of line strength. Insisting on 
ionization equilibrium between Fe I and Fe II allowed for the 
simultaneous determination of log g with T e ff and micro turbu- 
lence (£). Errors in T e s are ~ 50 K. We analyzed a reflected 
light spectrum of the asteroid Iris in order to obtain an instru- 
mental correction to the solar loge(Fe). Using this correction, 
we obtained [Fe/H] = 0.11 ±0.03 (s.e.m. 9 ) for the Praesepe 
cluster. 

2.5.2. Metallicity and Reddening Estimates in the Literature 

For other clusters in this study, we adopted or calculated av- 
erage metallicities from high-resolution spectroscopy as sum- 
marized in Table 5. For the studies by Boesgaard & Friel 
in Praesepe and the Pleiades, we have recomputed the clus- 
ter averages from individual [Fe/H] estimates after excluding 
cluster nonmembers. The only existing measurement of the 
abundance of NGC 2516 using high-resolution spectroscopy 
is from Terndrup et al. (2002) which was derived from only 
two stars. 

Table 6 lists reddening estimates from a variety of sources. 
For the Pleiades, we adopted a uniform reddening inferred 
from the H I hole measurement by Breger (1986). Although 
small parts of the cluster show higher reddening (e.g., Eggen 
1950; Mitchell & Johnson 1957), the stars in these regions 
constitute only a small fraction of our Pleiades sample, and 
most of these were tagged as outliers from our filtering algo- 
rithm (§ 4.1). We compared magnitudes and colors derived 

9 Standard error of the mean. 



Reference 


[Fe/H] 


s.e.m. 


Praesepe 


Boesgaard & Friel" 


+0.09 


0.03 


This paper (adopted value) 


+0.11 


0.03 


Pleiades 


Cayrel et al. (1988) 


+0.13 


0.07 


Boesgaard & Friel b 


+0.03 


0.02 


King et al. (2000) 


+0.06 


0.05 


Weighted mean 


+0.04 


0.02 






(0.03) c 


M67 


Garcia Lopez et al. (1988) 


+0.04 


0.04 


Hobbs&Thorburn(1991) 


-0.06 


0.03 


Friel & Boesgaard (1992) 


+0.02 


0.12 


Tautvaisiene et al. (2000) 


-0.03 


0.01 


Yong et al. (2005) 


+0.02 


0.02 


Randich et al. (2006) 


+0.03 


0.01 


Weighted mean 


+0.00 


0.01 






(0.03) c 


NGC 2516 


Terndrup et al. (2002) 


+0.01 


0.07 



a From Boesgaard (1989) and Friel & Boes- 
gaard (1992) after excluding known nonmem- 
bers. 

b From Boesgaard (1989) and Boesgaard & 
Friel (1990) after excluding known nonmem- 
bers. 

c Standard deviation of the measurements. 



from individual E(B-V) corrections for 157 stars (Breger 
1986;Stauffer&Hartmann 1987; Soderblomet al. 1993) with 
those from the uniform reddening, but the differences were 
negligible in our MS fitting range (§ 4.2). The reddening mea- 
surements for the Pleiades from other studies in Table 6 likely 
represent a cluster average including the CO cloud region. 

3. CALIBRATED ISOCHRONES 



Distances from Main Sequence Fitting. Ill 



7 



TABLE 6 

Cluster Reddening in the Literature 



Reference 


Method 


E(B-V) 


s.e.m. 


Praesepe 


Mermilliod (1981a) 


UBV 


0.00 


0.01" 


Nicolet (1981) 


Geneva 


0.011 


0.005 


Nissen (1988) 


uvby-0° 


0.007 


0.002 


Weighted mean 




0.007 


0.002 








(0.0()3) c 


Pleiades 


Mermilliod (1981a) 


UBV 


0.04 


O.OP 


Nicolet (1981) 


Geneva 


0.062 


0.005 


Breger (1986) d 


Sp-type, uvby-0° 


0.044 


0.003 


Nissen (1988) 


uvby-0° 


0.039 


0.005 


Weighted mean 




0.046 


0.005 








(0.0()9) c 


Breger (1986) [adopted value] e Sp-type, uvby-fS 


0.032 


0.003 


M67 


Taylor (1980) 


Multiple approaches 


0.046 


0.006 


Eggen(1981) 


uvby-0° 


0.050 


0.013 


Janes & Smith (1984) 


DDO 


0.056 


0.006 


Burstein etal. (1986) 


BH map f 


0.0358 


0.005 


Nissen et al. (1987) 


uvby-/3 b 


0.032 


0.006 


Hobbs&Thorburn(1991) 


T efi -(B-V) 


0.065 


0.018 


Montgomery et al. (1993) 


UBV 


0.05 


0.01 


Schlegel et al. (1998) 


Dust map 


0.032 


().005 h 


Weighted mean 




0.041 


0.004 








(0.010) c 


NGC 2516 


Eggen (1972) 


UBV 


0.125 


0.025 


Feinstein et al. (1973) 


Sp-type, UBV 


0.116 


0.004 


Snowden (1975) 


uvby-0 


0.120 


0.014 


Mermilliod (1981a) 


UBV 


0.11 


O.OP 


Nicolet (1981) 


Geneva 


0.114 


0.012 


Eggen (1983) 


uvby-0 3 


0.118 


0.004 


Verschoor & van Genderen (1983) VBLUW 


0.127 


0.005 


Nissen (1988) 


uvby-0° 


0.109 


0.007 


Dachs & Kabus (1989) 


Sp-type, UBV 


0.12 


0.004 


Sung et al. (2002) 


UBV 


0.112 


0.006 


Weighted mean 




0.117 


0.002 








(0.0()5) c 



"Assigned value in this paper. 

b E(B-V) = E(b -y)/0J4 (Cousins & Caldwell 1985). 
■■'Standard deviation of the measurements. 
d Average over the cluster field. 
e Average over the HI hole (see text). 
f Burstein & Heiles (1982). 

g Corrected for the systematic difference with Schlegel et al. (1998) of 0.02 mag. 
h Quoted precision (16%). 



We used the YREC to construct stellar evolutionary tracks 
at -0.3 < [Fe/H] < +0.2 in 0.1 dex increments. At each metal 
abundance, we ran a grid of masses from 0.2 to 8 M Q in 
0.05- IMq increments depending on the stellar mass. The 
scaled solar abundance mix from Grevesse & Sauval (1998) 
and the helium enrichment parameter AY / AZ =1.2 from 
non-diffusion models were used; see Sills et al. (2000) and 
Paper I for detailed information. The tracks were interpolated 
to generate theoretical isochrones at stellar ages from 20 Myr 
to 4 Gyr. 

Stellar luminosities and T e g were initially converted to V, 
B-V, V-Ic, and V-K s from the relations by Lejeune et al. 
(1997, 1998). To obtain finer grids of the isochrones, we in- 
creased the My spacing of the isochrone points using a cu- 
bic spline to match the AMy = 0.05 spacing of our tabulated 
color- r e ff corrections in Paper II. We then linearly interpolated 
the isochrones in both [Fe/H] and age at a fixed M v . We fi- 
nally applied the empirical color corrections, as defined in Pa- 
per II and updated in the Appendix of this paper. We assumed 



that the correction table is valid over the entire model grid in 
[Fe/H] and age. 
The isochrones constructed in this way are available online. 

10 



4. DISTANCES AT FIXED REDDENING 

We now proceed to derive the cluster distances, adopting 
the reddening values from Table 6. We do this in two ways, 
either holding the metallicity fixed at the spectroscopic values 
(Table 5) or determining a photometric metallicity that brings 
the distances from different CMDs into agreement. In either 
case, we show that the calibrated isochrones improve the in- 
ternal precision of distance estimates from the three CMDs 
with B-V, V-Ic, and V-K s as color indices, and V as a lumi- 
nosity index; hereafter (B-V,V), (V-I C ,V), and (V-K S ,V), 
respectively. 

10 See http://www.astronomy.ohio-state.edu/iso/. 
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4.1. Photometric Filtering 

Before fitting isochrones, we identified and removed stars 
that are far from the MS in comparison to the sizes of 
their photometric errors. Such stars could be either fore- 
ground/background objects or cluster binaries that stand off 
from the MS. 

The filtering process iteratively identifies the MS as the lo- 
cus of points of maximum density on the CMDs, indepen- 
dently of the isochrones. As an initial step, we identified 
the MS by hand, and then removed stars more than 1 mag 
away in V. This was necessary in particular for the analysis 
in NGC 2516, which shows many faint background stars in 
its CMDs (see below) that would complicate finding the den- 
sity maxima. The generous rejection criterion of 1 mag was 
chosen so that all cluster binaries would be preserved at this 
stage. 

At each step of the iteration, data points on each CMD were 
sorted by V magnitude into non-overlapping groups: each 
group contained \/N points, where N is the number of data 
points within a color range used in the filtering. The median 
colors in each group were computed, and then the resulting 
run of points (V vs. color) was made into a smooth curve by 
averaging each point with the linear interpolation of adjacent 
points. 

With the above trial MS, we computed x 2 from each CMD 

as 

N N (Ay,-) 2 

where AV, is the V difference between the i' h data point and 
MS at the color of the star and ay , and er c , are the photomet- 
ric errors in V and color, respectively. The error in the color 
contributes to the error in AV, by the slope of the curve, 7,, 
which was evaluated at the star's color. Because the above 
constructed curve would not precisely trace the observed MS 
in the presence of outliers, we added <jo in quadrature to the 
propagated photometric errors in the denominator. We ad- 
justed the value of ctq so that the total x 2 is equal to N. Ini- 
tially, we rejected all data points as the MS outliers if the xj 
(the individual contribution to x 2 ) is greater than 9 (corre- 
sponding to a 3 a outlier). We repeated adjusting <tq and re- 
jecting outliers with the reduced set of data points until there 
remained no point with xj greater than the above threshold 
value. We combined the results from all three CMDs and re- 
jected stars if they were tagged as an outlier from any one of 
the CMDs. 

We imposed an additional condition for the conver- 
gence of the algorithm because cluster binaries and fore- 
ground/background stars are often substantial populations 
near the MS. We compared the rms of AV, with its median 
of all absolute deviations (MAD). The MAD of a quantity x, 
is defined by 

MAD = 1 .48 3 median( |x, - median(x,) | ) , (2) 

where the correction factor 1.483 makes the estimator con- 
sistent with the standard deviation for a normal distribution 
(Rousseeuw 1990). The MAD is a more robust estimator 
of the dispersion than the standard deviation in the presence 
of outliers. Therefore, we reduced the size of the x} thresh- 
old value if the fractional difference between the rms and the 
MAD of AV, is larger than 5%, assuming median(AV,) = 0. 
We repeated the above filtering steps by constructing a new 
MS curve from the remaining data set. 




0.0 0.5 1.0 1.5 
B - V 

FIG. 4. — Result of the filtering process in Praesepe. All known cluster 
binaries were included to test the filtering algorithm. The plus signs are re- 
jected points from the filtering, and open circles are those remaining. Filled 
circles are known binaries that remained after the filtering. Stars brighter than 
V = 7 were excluded before the filtering algorithm was applied. 



Figure 4 illustrates the result of the filtering algorithm in 
Praesepe. Here we include all known cluster binaries but ex- 
clude stars with a low membership probability (§ 2.4). We 
used all three CMDs in the filtering, but only the (B - V, V) 
diagram is shown. The plus signs are rejected points from 
the filtering, and open circles are those remaining. Some of 
the stars on the MS were rejected because they were filtered 
in other CMDs. Although many binaries remain after the fil- 
tering (filled circles), their proximity to the MS would have 
a minor impact on the derived distance. The bias in the dis- 
tance due to the remaining binaries is discussed in § 6.2 from 
artificial cluster CMD tests. 

4.2. Isochrone Fitting and the Photometric Metallicity 

The cluster distances were found by fitting isochrones over 
a range of metallicity with the adopted age of each cluster 
(§ 2.1). The isochrones were reddened to the E(B-V) values 
in Table 6. We adopted a reddening prescription for the broad- 
band colors by Bessell et al. (1998), who have assumed the ex- 
tinction law of Mathis (1990). Their reddening formulae were 
computed from ATLAS9 synthetic stellar photometry for a 
wide range of T e ff, and include color-dependent reddening re- 
sulting from shifts in the effective wavelengths of broadband 
filters. The Bessell et al. formulae give reddening and ex- 
tinction values for E(B-V) = 0.30, so we linearly rescaled 
them according to the assumed cluster reddening. The color 
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FIG. 5. — Example of determining a photometric metallicity in the Hyades. 
Open circles are the distances derived from CMDs with B — V, V — Iq, and 
V -K s as a function of isochrone metallicity. The solid lines connecting these 
points are least-squares fits, and the labels to the left of the lines indicate the 
corresponding color index. An age of 550 Myr and E(B-V) = 0.000 were 
assumed. The photometric metallicity is defined as the average [M/H] where 
the B — V line crosses the V—Ic andV-K s lines. 

transformation by Carpenter (2001) 11 was used to compute 
the reddening in V — K s . For zero-color stars, we found 
R v = A V /E(B-V) = 3.26, R V i = E(V -I C )/E(B-V) = 1.32, 
and R VK = E(V - K s ) / E(B -V) = 2.91. For stars in the middle 
of our MS-fitting range, (B-V)o = 0.8, we found R v = 3.44, 
R V i= 1.37, and R VK = 3.04. 

For a given isochrone, we computed each individual star's 
distance modulus (/!,), and defined the cluster distance mod- 
ulus in each CMD as the "unweighted" median of /i,, i.e., 
(m-M)o = median(^,). We computed the fitting error in the 
distance modulus on each CMD as 



C(m-M) - max 



MADQm) 



(3) 



phot '' Vn 

where N is the total number of data points used in the fit. Here 



1 



-phot 



i 



(4) 



where the same notation is used as in equation (1) except that 
7, is the isochrone slope at the color of the star. We further 
computed a "weighted" mean and a "weighted" median of 
the distance modulus, and discuss the difference between the 
three distance estimates along with other systematic errors in 
§6.1. 

For a given set of [M/H], age, and E(B - V), we fit 
isochrones in (B-V,V), (V - I C ,V), and (V - K S ,V) over 
0.4 < (B-V)o < 1.3. This color range is where the Hyades 
calibration is most reliable (Paper II). There were only a 
few Hyades members blueward of the range, and the cool 
end was set by the magnitude limit of the Hipparcos mis- 
sion. The range corresponds to 0.48 < (V—Ic)o < 1 .48 and 
0.98 < (V-K s ) < 3.16. For non-zero reddening, the fitting 
ranges were made correspondingly redder. 

We define the photometric metallicity [M/H] £ as the one 
that brings distances from three CMDs into statistical agree- 
ment (e.g., Pinsonneault et al. 1998; Stello & Nissen 2001; 

11 Updated color transformations 2MASS All-Sky Data Release can be 
found at 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4b.html. 



Terndrup et al. 2002). Figure 5 shows how this process works 
for the Hyades, our calibrating cluster. The open circles dis- 
play the derived distances as a function of isochrone metal- 
licity in each of the three CMDs with E(B-V) = 0.000 at 
an age of 550 Myr. The lines connecting these points are 
least-squares fits, and the labels to the left of the lines indi- 
cate the corresponding color index. The slope in the (B—V,V) 
is larger than in the other two CMDs, indicating a greater 
sensitivity of the isochrone luminosity to the metal abun- 
dance. We define [M/H] £ as the weighted average of the two 
metallicities at which the (B - V, V) distance agrees with that 
from (V-Ic,V) and from (V-K S ,V), respectively. For the 
Hyades in Figure 5, we derive the photometric metallicity of 
[M/H] £ = +0.13 ± 0.02, which is naturally the same as the 
originally assumed value in the isochrone calibration. 

With this definition of the photometric metallicity, we can 
derive two distances at the adopted reddening. The first of 
these, designated as (m-M)o^, is the weighted average dis- 
tance modulus from the three CMDs at the spectroscopic 
[Fe/H]. The second distance, (m-M)o ; „, is the value deter- 
mined at the photometric metallicity, [M/H] £ . 

4.3. Results 

Figures 6-9 display the CMDs for each cluster. For Prae- 
sepe and the Pleiades, known binaries and nonmembers were 
excluded before we applied the photometric filtering. Stars 
that remained after the filtering are shown as open points, and 
the rejected stars are shown as plus signs. With the exten- 
sive binary and membership information, the CMDs of Prae- 
sepe (Fig. 6) are dominated by single cluster members, and 
there are few stars rejected by the photometric filtering. The 
CMDs of the Pleiades (Fig. 7) show that the filtering routine 
identified likely cluster binaries effectively. The single cluster 
members from S04 are shown in the CMDs of M67 (Fig. 8), 
none of which were rejected by the filtering algorithm. The 
CMDs of NGC 25 16 (Fig. 9) are from JTH01 and have a large 
number of foreground/background stars as expected in an area 
survey (for clarity, the point size for the rejected stars was re- 
duced on the CMDs as compared to the other plots). 

Isochrones at the spectroscopic metallicities are overlaid as 
solid lines in each figure, and in most cases they are in excel- 
lent agreement with the observed shapes of the MS. In partic- 
ular, the match to the Praesepe MS (Fig. 6), which has an age 
and metallicity nearly identical to the Hyades, is good over 
almost 7 mag in V. We interpret this agreement as indicating 
that the Hyades-based corrections to the color- r e ff relations 
defined in Paper II apply to all clusters, or at least to those 
with metallicities not too different from that of the Hyades. 

There are a few cases in which the match to the MS shape is 
not as good as it is in Praesepe. In the Pleiades (Fig. 7), stars 
with B-V > 0.9 are considerably bluer than the isochrone in 
(B-V,V), although not in the other two CMDs. This "blue K 
dwarf" phenomenon was discussed by Jones (1972), Landolt 
(1979), and van Leeuwen et al. (1987), and at greater length 
by Stauffer et al. (2003), who attributed it to stellar temper- 
ature inhomogeneities caused by cool spots and plage areas 
in rapidly rotating young stars. From the JTH01 photometry, 
it would appear that NGC 2516, which is about 40% older 
than the Pleiades, does not share this phenomenon with the 
Pleiades (see Fig. 9). However, we show in § 7.1 that the 
usage of the S02 data comes to a different conclusion with 
supporting evidence from stellar rotational velocities. The 
data in (V-/ C ,V) for all clusters are somewhat bluer than 
the isochrone redward of V-Ic ~ 1-6; this indicates a pos- 
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FIG. 6. — CMDs of Praesepe, after excluding known binaries and stars with low membership probability. Plus signs are photometrically rejected data points 
from the filtering algorithm, and open circles are those remaining. The solid lines are empirically calibrated isochrones with spectroscopic metallicity (Table 5), 
which were adjusted for the reddening with the literature value (Table 6). Fitting ranges are shown as horizontal bars. The arrow denotes the direction of 
reddening vectors. 




FIG. 7. — Same as Fig. 6, but for the Pleiades. As discussed by Stauffer et al. (2003), the Pleiades K dwarfs are bluer than the given isochrone in (B-V,V), but 
not in the other two CMDs (see text). The distance modulus in (B - V, V) was derived at 0.4 <(B- V) < 0.8 as shown by the shorter horizontal bar. 
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FIG. 9. — Same as Fig. 6, but for NGC 2516 with photometry from JTH01. The fitting range in (B-V, V) was truncated in the red part as in the Pleiades CMD 
(see text). 
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FIG. 10. — Distance vs. isochrone metallicity for Praesepe with the Leje- 
une et al. (1997, 1998) color-r e ff relation (top) and with our new empirical 
correction (bottom). The reddening value in the literature (Table 6) was as- 
sumed. Circles and lines have the same meaning as in Fig. 5. Vertical lines 
show the spectroscopic metallicity and its 1 <r error (Table 5). 



sible systematic error in the Hyades calibration for these red 
stars, probably resulting from a paucity of Hyades stars with 
accurately measured parallaxes in this color range. In M67 
(Fig. 8), there may be a larger deviation because the color cal- 
ibration employed by S04 did not extend as far to the red as 
this. 

The distance was independently derived from each CMD 
over the color range shown as a horizontal bar in each panel. 
The blue end of the color ranges for M67 was increased to 
(B-V)o = 0.7 to avoid the steeply rising parts of the MS, as the 
resulting distance would be sensitive to the choice of cluster 
age. For the Pleiades and NGC 2516, the red end of the color 
range in (B-V,V) was decreased to (B-V)q = 0.8 to avoid 
systematic errors from the blue K dwarf phenomenon. 

In Table 7 we assemble our derived cluster distances and the 
photometric metallicities. The first column lists the source of 
photometry, including particular selections of subsamples of 
the data. We applied the photometric filtering to each sub- 
sample before fitting isochrones. The second through fourth 
columns display the distance modulus found from each CMD, 
where the spectroscopic metallicity and reddening were taken 
from Tables 5 and 6, respectively. The weighted average dis- 
tance, (m-M)o.s, is shown in the fifth column. The last two 
columns display the photometric metallicity, [M/H] £ , and the 
distance modulus, (m-M)o,£, at this metallicity. The error in 
(m-M)o i j is either the propagated one from the errors of the 
individual distances or the s.e.m. of these distance estimates, 
whichever is larger. The error in (m-M)o.£ additionally in- 
cludes the propagated error from [M/H] £ . 
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FIG. 11.— Same as Fig. 10, but for M67 with MMJ93 (top) and S04 
(bottom) photometry. The stars designated as single cluster members by S04 
were used in both cases. 

The standard errors of the average distance modulus at the 
spectroscopic metallicities are typically on the order of 0.02 
mag (i.e., the individual distances in the three CMDs are con- 
sistent to <~ 2%). In Figure 10 we show that this consistency 
is primarily the result of the empirical color- T e ff corrections 
derived in Paper II. The top panel shows the derived distances 
to Praesepe as a function of the metallicity for isochrones 
that do not incorporate these corrections, while the bottom 
panel shows what happens when the corrected isochrones are 
employed. The vertical lines show the metallicity estimate 
for Praesepe with its 1 o error (§ 2.5.1). The internal preci- 
sion of the distance estimation with isochrones employing the 
Hyades empirical calibration is about a factor of 5 better than 
with the uncorrected isochrones. 

The distances and photometric metallicities in Table 7 are 
generally not sensitive to the selection of subsamples in each 
cluster. For example in Praesepe and the Pleiades, leaving 
known binaries in the sample before filtering changed the 
value of (m-M)o.s by only about 0.01 mag, confirming that 
the filtering algorithm removed most of the binaries that are 
brighter than the MS (e.g., Fig. 4). In M67, we compared the 
distances from the MMJ93 data, but also using their photom- 
etry for only those stars identified as single members by S04; 
we found a negligible difference. In NGC 25 16, we computed 
the distances from the full JTH01 or S02 catalogs, or select- 
ing only those stars identified as radial velocity (RV) members 
(Terndrup et al. 2002) or X-ray detected sources (Damiani 
et al. 2003); again making that selection reduced the distance 
modulus by only about 0.01 mag. 

Figures 11-13 display the derived distances as a function 
of metallicity for M67, NGC 2516, and the Pleiades, respec- 
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TABLE 7 

MS Fitting Distance and Photometric Metallicity at Literature E(B-V) 







(m-M)o 






[M/H] £ 




Data / Constraint 


U»-V,V) 


(V-lc.V) 


(V-K S .V) 


(m-M) ,s 


(m-M) ,£ 


Praesepe 


known binaries excluded 
known binaries included 


6.267 ±0.008 
6.262 ± 0.008 


6.291 ±0.008 
6.279 ± 0.008 


6.336 ±0.007 
6.327 ± 0.007 


6.303 ±0.021 
6.293 ±0.020 


+0.207 ±0.030 
+0.186 ±0.030 


6.395 ±0.030 
6.367 ±0.031 


Pleiades 


known binaries excluded 
known binaries included 


5.658 ±0.010 
5.660 ±0.009 


5.623 ±0.013 
5.639 ±0.013 


5.669 ±0.017 
5.677 ±0.016 


5.649 ±0.013 
5.657 ±0.009 


+0.021 ±0.032 
+0.035 ±0.028 


5.627 ±0.039 
5.651 ±0.036 


M67 


S04 

MMJ93 

MMJ93 / S04 stars 2 


9.631 ±0.008 
9.630 ±0.010 
9.612 ±0.013 


9.601 ±0.007 
9.505 ±0.009 
9.529 ±0.010 


9.654 ±0.011 
9.626 ±0.008 
9.631 ±0.011 


9.622 ±0.015 
9.584 ±0.042 
9.584 ±0.033 


-0.002 ± 0.062 
-0.032 ±0.183 
+0.017±0.138 


9.620 ± 0.054 
9.562±0.163 
9.599 ±0.109 


NGC 2516 


JTHOl 

JTHOl /RV + X-ray b 
S02 

S02 / RV + X-ray b 


8. 125 ±0.007 
8.115±0.010 
8.098 ±0.014 
8.081 ±0.013 


8.040 ± 0.009 

8.041 ±0.014 
8.035 ±0.013 
7.960 ±0.022 


8.100 ±0.007 
8.077 ±0.021 
8.064 ±0.012 
8.063 ±0.022 


8.096 ±0.023 
8.088 ±0.023 
8.065 ±0.017 
8.052 ±0.034 


-0.052 ± 0.053 
-0.074 ±0.030 
-0.042 ± 0.024 
-0.059 ±0.081 


8.034 ±0.058 
7.995 ±0.037 
8.015 ±0.026 
7.979 ±0.092 



NOTE. — All subsamples of data were filtered independently before fitting isochrones. 
"Single cluster members from S04 with MMJ93 photometry. 

b RV members from Terndrup et al. (2002), and X-ray detected sources from Damiani et al. (2003). 
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FIG. 12.— Same as Fig. 10, but for NGC 2516 with JTHOl (top) and 
S02 (bottom) photometry. The fitting range in (B-V.V) was restricted to 
0.4 < (B-V)o < 0.8 as in the case of the Pleiades. The dashed line shows 
the result when the full fitting range was used in this color. 

tively. The vertical lines on each plot show the average [Fe/H] 
and its 1 a error from Table 5. The distances for M67 (Fig. 11) 
are significantly less consistent with the MMJ93 photometry 
than from the S04 data. This is also seen in the larger errors 
in [M/H] £ and (m-M)o,£ in Table 7. This is mainly because 
the V-Ic colors in MMJ93 are (on average) 0.022 mag redder 
than in S04, which will reduce the distance in the (V - Ic, V) 
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FIG. 13. — Same as Fig. 10, but for the Pleiades. Horizontal lines represent 
the average geometric distance and its 1 a error from Table 1. 

by about 0.07 mag since the average slope of the MS in this 
color is about 3. 

In Figure 12 we explore the effects of the blue K dwarfs 
in NGC 2516 by comparing the two photometric data sets. 
The distances in (B-V,V) derived over the truncated color 
range are shown as a solid line, while the dashed line displays 
the results for fits over 0.4 < (B-V) < 1.3. The fits in the 
other two CMDs were done over the full color ranges. When 
the truncated color range is used, both JTHOl and S02 data 
give distances and photometric metallicities that are mutually 
consistent. However, fitting on the S02 data over the full range 
of colors in B-V reduces the photometric metallicity by <~ 
0.12 dex and the average distance at this metallicity by ~ 0.09 
mag, while the results from the JTHOl data remain consistent 
within 1 a independently of the fitting range. We discuss the 
blue K dwarfs phenomenon in NGC 2516 at greater length in 
§7.1. 

In Figure 13 we display the fits for the Pleiades, and also 
plot the average geometric distance and its 1 a error from 
Table 1 as the horizontal lines. The agreement between the 
MS-fitting and geometric distances is superb. Our distance 
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for Praesepe at the spectroscopic metallicity also agrees with 
existing geometric distance determinations: Gatewood & de 
Jonge (1994) found (m-M) = 6.42 ±0.33 from ground- 
based parallaxes, and Loktin (2000) obtained 6.16 ±0.19 
from the moving cluster method. In addition, there are two 
determinations of the distance to Praesepe from Hipparcos: 
(m—M)o = 6.37 ± 0. 15 (van Leeuwen 1999) and \ 3 2 (Ro- 
bichon et al. 1999). We compare our derived distances with 
the Hipparcos measurements in § 7.2. 

5. SIMULTANEOUS DETERMINATION OF CLUSTER PARAMETERS 

In the previous section, we showed that the cluster metal- 
licities can be obtained from the MS fitting because the vari- 
ous color indices have different sensitivities to the metallicity. 
However, we can further constrain the reddening because the 
slope of the MS becomes shallower below My ~ 5.5, espe- 
cially in (V-Ic,V) or (V-K S ,V). Furthermore, the reddening 
vectors are not parallel to the MS as shown in Figure 6, so 
the derived distances depend on the reddening with different 
degrees of sensitivity in the three CMDs. 

To determine the metallicity, reddening, and distance simul- 
taneously, we minimized 



Xtot : 



3 Nj 

EE 

;■=! ;=1 



(5) 



where the same notation is used as in equation (1) except that 
the subscript j is used to indicate each of the three CMD. The 
Nj is the number of data points used in the fit on the /* CMD, 
and fiij are the individual distance moduli of the i th point. The 
Ji is the weighted average distance modulus from the three 
CMDs. 

The quantity fj in the denominator of equation (5) is a 
renormalization factor to take into account underestimation 
of the photometric errors. The value of fj was determined 
by comparing the residuals about the best-fitting isochrone at 
the minimum x 2 0t (see below) to the propagated distance error 
from the photometric errors (eq. [4], er p h t)- This yields 



fj = max 



MAD(p u ) 1 



Nj crp hot j ' 



1 



(6) 



This renormalization was required in particular for the S04 
data, which had many more repeat measurements in V and 
I c than in B, so the quoted errors are much smaller in V - Ic 
than in the other colors. For these data, we found / w 6 in 
(V-Ic,V), indicating that the errors are greatly underesti- 
mated at least when compared to the scatter about the best- 
fitting isochrone. Without the renormalization factor, the ma- 
jority of the weight in the x 2 t would be given to the (V-I c , V) 
CMD, and the process would almost entirely ignore the infor- 
mation on the other CMDs. 

We searched for the minimum Xtot m tne plane of [M/H] 
versus E(B-V) using a downhill simplex method (Press et al. 
1992). The average distance modulus in equation (5) was 
determined at each [M/H] and E(B-V). Because the color 
ranges used in the fitting depend on the adopted E(B-V), we 
solved for the minimum x 2 t at tne E(B-V) derived from the 
previous iteration to keep the same number of data points for 
each iteration. Since we defined Xm w i tn respect to the av- 
erage distance modulus from the three CMDs, the minimum 
X 2 ot yields a set of [M/H] and E(B-V) that best describes all 
of the CMDs simultaneously. 



Figure 14 shows likelihood contours in the [M/H] versus 
E(B-V) plane for all clusters in this paper and for the Hyades. 
The x 2 surfaces were smoothed with a Gaussian kernel. Con- 
tours are shown at Ax^t = 3.53, 8.02, and 14.2 relative to the 
minimum value of %f ot , which correspond to 68.3%, 95.4%, 
and 99.7% confidence levels (1, 2, 3 a) for the 3 degrees of 
freedom, and the errors in both quantities were assumed to 
be normally distributed without correlation. The vertical lines 
on each panel display the spectroscopic metallicity and its 1 
a error, while the horizontal lines show the average reddening 
in the literature with its 1 a error. The error bars in each panel 
show the size of systematic errors, which is discussed in the 
next section. 

The solution for Praesepe (Fig. 14, middle left panel) yields 
a value for the reddening that is slightly negative and a metal- 
licity that is <~ 2 a higher than found spectroscopically. The 
negative reddening may be a consequence of zero-point er- 
rors in the photometry or else from an error in the Hyades 
reddening, which we assumed in Paper II as identically zero. 
The study of Taylor (1980), for example, indicates E(B-V) = 
+0.003 ± 0.002 for the Hyades. So, if the Praesepe redden- 
ing was slightly less than the non-zero Hyades reddening, the 
MS-fitting method would yield a negative value. 

In Figure 14 two solutions are shown for M67 (top mid- 
dle panel): the contours with the higher E(B-V) are from 
the MMJ93 data, and the others are from the S04. The large 
difference in E(B-V) obtained from these two samples is pri- 
marily due to the difference in the V— Ic values (Table 2). We 
illustrate this further by showing the x 2 contours that result 
when only the B - V and V - K s data are included. The mid- 
dle panel shows the likelihood contours for the MMJ93 data, 
while the bottom middle panel is for the S04. Here the two 
studies yield the same values of metallicity and reddening, al- 
though with larger errors. We also display two solutions for 
NGC 2516 in the bottom left panel of Figure 14. The con- 
tours with the higher reddening are from the SBLL03, and the 
others are from JTH01 photometry. 

The top right panel of Figure 14 shows the solution for the 
Pleiades, while the panel below that shows the solution for 
the metallicity and reddening that result if the average geo- 
metric distance from Table 1 is taken as a prior. Contours in 
this panel are 1,2, and 3 a confidence levels for the 2 degrees 
of freedom. This demonstrates that good photometry and par- 
allaxes together can provide strong constraints on the metal 
abundance of a cluster. 

Table 8 lists the results from the % 2 solutions for vari- 
ous combinations of the data listed in the first column. The 
[M/H] x 2 and E(B-V) x i denote our derived values of metal- 
licity and reddening, respectively, and are shown in the third 
and fourth columns. Errors in these quantities are the sizes of 
the 1 (7 contours in Figure 14. The (m-M) Q x 2 in the second 
column is the average distance from the three CMDs. The 
errors in metallicity and reddening were propagated into the 
error in distance. We compare these values with previous es- 
timates in the literature after a discussion of systematic errors 
in the next section. 

6. SYSTEMATIC ERRORS AND THE ACCURACY OF MS FITTING 

In the previous two sections we showed that MS fitting can 
be performed with high internal precision, resulting in errors 
in distance moduli of 0.02 mag (i.e., 1% in distance) at the 
spectroscopic metallicity (Table 7). When spectroscopy is not 
available, photometry in BVIcK s alone can be used to derive 
distances with precision of 1 % - 3%, metallicities with 0.02 - 
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FIG. 14. — Likelihood contours in [M/H] vs. E(B-V) shown at Axj? ot = 3.53, 8.02, and 14.2 (1, 2, and 3 cr for 3 degrees of freedom) relative to the minimum 
value of x^ ot (eq. [5]). Those for 2 degrees of freedom are shown for the Pleiades at the geometric distance (middle right). For NGC 2516 (bottom left), contours 
with the higher E(B-V) are from the S02, and the others are from the JTH01 photometry. For M67 (top middle), contours with the higher E(B-V) are from 
the MMJ93, and the others are from the S04 photometry. Also shown are contours for the MMJ93 (middle) and the S04 (bottom middle) in the case where only 
(B-V,V) and (V-K S ,V) are used. Error bars represent 1 o systematic errors listed in Tables 9-12 excluding fitting errors. The spectroscopic metallicities and 
reddening values in the literature with 1 o errors (Tables 5, 6) are shown with vertical and horizontal lines, respectively. 



0.04 dex error, and reddening estimates with better than 0.01 
mag error (Table 8). However, various systematic errors in 
the MS fitting should be considered in the final error budget 
such as the photometric zero-point errors or biases introduced 
by the cluster binaries and foreground/background stars. In 
this section we estimate these systematic errors and show that 
they often exceed the internal precision of the MS fitting by 
factors of 2 - 3. 

6.1. Errors from Input Quantities 

In Tables 9-12 we list the sources of several systematic er- 
rors and their contributions to errors in the distance, photo- 
metric metallicity, and reddening estimates for each cluster. 
The first column displays the source of errors, and the second 
column shows the size of the errors adopted for each cluster. 
The third column lists error contributions to (m-M)o.s when 
both metallicity and reddening are held fixed. The fourth and 
fifth columns contain errors in (m-M)o,£ and [M/H] £ when 
the photometric metallicity is determined at a fixed reddening. 
The sixth through eighth columns list errors in (m-M) Q x i, 



[M/H] x 2, and E(B-V) x i from the x 2 minimization. At the 
bottom of the table, we list the "total error," computed as the 
quadrature sum of all of the error contributions. These are the 
errors that are shown as error bars along with the x 2 contours 
in Figure 14. The fitting errors (i.e., the internal precision of 
MS fitting) are not included in these sums but instead are sep- 
arately tabulated below. 

In the first column of these tables, the line marked [Fe/H] 
denotes the error in the adopted cluster metallicity (Table 5). 
This affects the error in (m-M)o.s only. The errors in the 
adopted reddening (Table 6) are shown in the line marked 
E(B-V). These errors are propagated into the determina- 
tion of the photometric metallicity [M/H]e and the distance 
(m-M)o.E, but not the values determined from the x 2 mini- 
mization. The remaining errors will contribute to the system- 
atic errors of all of the parameters found via MS fitting. 

We chose an error of 30% in the adopted age of each 
cluster from the consideration of previous age estimates for 
the Pleiades. The cluster age is about 100 Myr from the 
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TABLE 8 
X 2 Solutions 



Data / Constraint 


(m-M) 0x 2 


[M/H] x2 


E(B-V) X 2 


Praesepe 


known binaries excluded 
known binaries included 


6.365 ±0.025 
6.357 ±0.026 


+0.202 ±0.020 
+0.192 ±0.020 


-0.003 ±0.005 
0.002 ±0.004 


Pleiades 


known binaries excluded 
At geometric distance a 
known binaries included 
At geometric distance c 


5.693 ±0.053 
5.632 ±0.017 b 
5.667 ±0.054 
5.632 ±0.017 b 


+0.079 ±0.040 
+0.058 ±0.012 
+0.059 ±0.040 
+0.049 ±0.012 


0.032 ±0.007 
0.024 ±0.004 
0.027 ±0.007 
0.025 ±0.004 


M67 


S04 

MMJ93 

MMJ93 / S04 stars d 


9.558 ±0.034 
9.603 ±0.051 
9.638 ±0.053 


-0.043 ±0.032 
-0.061 ±0.035 
+0.019 ±0.050 


0.024 ±0.006 
0.078 ±0.006 
0.061 ±0.010 


NGC 2516 


JTH01 

JTH01/RV + X-ray c 
S02 

S02 / RV + X-ray 6 


8.001 ±0.038 
7.939 ±0.059 
8.064 ±0.038 
8.038 ±0.091 


-0.076 ±0.020 
-0.133 ±0.040 
-0.055 ±0.030 
-0.050 ±0.060 


0.1 12 ±0.007 
0.120 ±0.010 
0.1 39 ±0.007 
0.134±0.015 



NOTE. — All subsamples of data were filtered independently before fitting 
isochrones. 



"Solution at the average geometric distance after excluding known binaries. 
b Average geometric distance from Table 1. 

c Solution at the average geometric distance with known binaries included in the 
data. 

''Single cluster members from S04 with MMJ93 photometry. 
e RV members from Terndrup et al. (2002), and X-ray detected sources from Dami- 
ani et al. (2003). 



TABLE 9 

MS Fitting Error Budget for Praesepe 



Source of Error 


AQuantity 


Adopted [Fe/H] 
A(m-M)o. s 


Adopted E(B-V) 


x 2 


Minimization 




A(m-M) ,£ 


LM/Hj £ 


A(m-M) 0x 2 


[M/Hj x2 


E(B-V) X 2 


[Fe/H] 


±0.03 


±0.029 












E(B-V) 


±0.002 


±0.005 


±0.013 


±0.007 








Age 


±30% 


±0.007 


±0.010 


±0.000 


=F0.014 


=F0.013 


q=0.004 


Helium (Y) . . . 


±0.009 


=F0.027 


T0.027 


±0.000 


q=0.027 


±0.000 


±0.000 


Calibration . . . 




±0.010 


±0.020 


±0.020 


±0.010 


±0.010 


±0.002 


Fitting method 




±0.006 


±0.014 


±0.013 


±0.002 


±0.004 


±0.001 


Rv 


±0.3 


±0.000 


±0.006 


±0.006 


±0.000 


=F0.001 


±0.000 


A'; ; 


±0.07 


±0.001 


±0.004 


±0.003 


T0.002 


T0.003 


±0.001 


Rvk 


±0.12 


±0.001 


±0.002 


±0.001 


±0.000 


±0.000 


±0.000 


AV 


±0.006 


±0.001 


=F0.011 


=F0.011 


T0.003 


T0.008 


±0.002 


A A. 


±0.007 


±0.006 


±0.019 


±0.013 


±0.010 


±0.010 


T0.002 


A(B-V) 


±0.004 


T0.006 


±0.038 


±0.044 


±0.014 


±0.027 


T0.002 


A(V-I) C 


±0.006 


T0.008 


T0.027 


T0.021 


T0.007 


T0.006 


±0.002 


Total 




±0.044 


±0.065 


±0.058 


±0.037 


±0.035 


±0.006 


Fitting 




±0.021 


±0.030 


±0.030 


±0.025 


±0.020 


±0.005 



MS turnoff using isochrones that incorporate convective over- 
shooting (Meynet et al. 1993), while ages from the lithium 
depletion boundary are about 125 Myr (Stauffer et al. 1998; 
Burke et al. 2004). However, distances, photometric metallic - 
ities, and reddening values are insensitive to age because the 
color range we chose for the fitting avoids the upper MS and 
the pre-MS (for clusters at least as old as the Pleiades). 

The helium abundance (Y) sensitively affects the isochrone 
luminosity (AM V /AY rs 3 at fixed T e ff). A shorter distance 
is obtained when a higher helium abundance is assumed. We 
adopted an error in Y of 0.009 for each cluster. This is a 1 a 
scatter in AY / AZ from the primordial helium abundance, the 
solar value, and that of the Hyades as inferred from the MS 
luminosity at the Hipparcos distance (see Paper I). The photo- 
metric metallicities and reddening values are not affected by 



the helium abundance since the error in Y equally changes the 
distances in the three CMDs. 

The calibration errors in the tables reflect errors in our 
adopted parameters of the Hyades: [Fe/H] = +0.13 ±0.01 
(Paulson et al. 2003) and the center-of-mass distance modulus 
of (m-M) = 3.33 ±0.01 (Perrymanet al. 1998). We adopted 
a 1 a error in reddening from Taylor (1980), E(B-V) = 
+0.003 ± 0.002, but assumed a zero reddening toward the 
Hyades (Paper II). The errors in [M/H] £ and (m-M)o,E fur- 
ther include the errors arising from the scatter about the linear 
relation between the two quantities (see Fig. 5). The error 
from the fitting method comes from experiments in which we 
used the weighted mean or weighted median distance modu- 
lus instead of the unweighted median value. 

The next three rows in the tables show the effects of er- 
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TABLE 10 

MS Fitting Error Budget for the Pleiades 



Source of Error 


AQuantity 


Adopted [Fe/H] 
Aim— M)a i 


Adopted E(B-V) 


x 2 


Minimization 




A(m-M) .E 


|M/H| £ 


A(m-M) n 2 


[M/HJ X 2 


E(B-V) 2 


[Fe/H] 


±0.02 


±0.023 












E(B-V) 


±0.003 


±0.007 


ztU.UIU 










Age 


±30% 


+0.006 


+0.020 


+0.011 


+0.004 


+0.002 


+0.001 


Helium (Y) . . . 


±0.009 


=F0.027 


T0.027 


±0.000 


=F0.027 


±0.000 


±0.000 


Calibration . . . 




±0.010 


±0.020 


±0.020 


±0.010 


±0.010 


±0.002 


Fitting method 




±0.006 


±0.038 


±0.033 


±0.027 


±0.014 


±0.005 


Rv 


±0.3 


T0.003 


±0.018 


±0.018 


±0.013 


±0.016 


T0.001 


Rvi 


±0.07 


±0.003 


±0.012 


±0.008 


±0.010 


±0.008 


to.ooi 


Rvk 


±0.12 


±0.002 


±0.008 


±0.005 


±0.007 


±0.005 


±0.000 


AV 


±0.012 


±0.008 


T0.009 


=F0.014 


T0.013 


=F0.018 


±0.000 


AK, 


±0.007 


±0.003 


±0.012 


±0.008 


±0.014 


±0.010 


±0.000 


A(B-V) 


±0.008 


T0.023 


±0.058 


±0.069 


±0.048 


±0.066 


T0.002 


A(V-l)c 


±0.009 


T0.013 


T0.037 


^0.021 


T0.034 


T0.030 


±0.006 


Total 




±0.048 


±0.092 


±0.087 


±0.076 


±0.080 


±0.009 


Fitting 




±0.013 


±0.039 


±0.032 


±0.053 


±0.040 


±0.007 



TABLE 11 
MS Fitting Error Budget for M67 



Adopted [Fe/H] Adopted E{B - V) \ Minimization 



Source of Error 


AQuantity 


A(m-M) 0]S 


A(m-M) 0£ 


[M/HJ £ 


A(m-M) 0x 2 


[M/Hj x 2 


U(B-V) x2 


[Fe/H] 


±0.01 


±0.008 












E(B-V) 


±0.004 


±0.008 


±0.016 


±0.009 








Age 


±30% 


±0.006 


±0.009 


±0.005 


±0.001 


T0.002 


=F0.002 


Helium (Y) . . . 


±0.009 


T0.027 


T0.027 


±0.000 


T0.027 


±0.000 


±0.000 


Calibration . . . 




±0.010 


±0.020 


±0.020 


±0.010 


±0.010 


±0.002 


Fitting method 




±0.004 


±0.011 


±0.016 


±0.001 


±0.001 


±0.000 


Rv 


±0.3 


T0.004 


±0.019 


±0.026 


±0.000 


±0.008 


T0.002 


Rvi 


±0.07 


±0.006 


±0.014 


±0.010 


±0.004 


±0.002 


to.ooi 


Rvk 


±0.12 


±0.001 


±0.007 


±0.007 


±0.002 


±0.001 


±0.000 


Ar 


±0.009 


±0.006 


T0.006 


T0.014 


T0.007 


T0.016 


±0.000 


A A. 


±0.007 


±0.002 


±0.011 


±0.011 


±0.012 


±0.012 


±0.000 


A(B-V) 


±0.005 


=F0.006 


±0.038 


±0.049 


±0.026 


±0.043 


T0.003 


A(V-/)c 


±0.011 


T0.022 


T0.045 


T0.023 


T0.003 


T0.009 


±0.014 


Total 




±0.040 


±0.076 


±0.070 


±0.041 


±0.050 


±0.014 


Fitting 




±0.015 


±0.054 


±0.062 


±0.034 


±0.032 


±0.006 


NOTE. — Estimated from S04 photometry. 



TABLE 12 

MS Fitting Error Budget for NGC 25 16 



Source of Error 


AQuantity 


Adopted [Fe/H] 
A(m-M)o. s 


Adopted E(B-V) 


x 2 


Minimization 




A(m-M) ,£ 


LM/HJ £ 


A(m-M) 0x 2 


[M/HJ X 2 


E(B-V) X 2 


[Fe/H] 


±0.07 


±0.072 












E(B-V) 


±0.002 


±0.006 


±0.005 


±0.000 








Age 


±30% 


±0.005 


±0.014 


±0.008 


±0.014 


±0.006 


±0.000 


Helium (Y) . . . 


±0.009 


=F0.027 


T0.027 


±0.000 


T0.027 


±0.000 


±0.000 


Calibration . . . 




±0.010 


±0.020 


±0.020 


±0.010 


±0.010 


±0.002 


Fitting method 




±0.009 


±0.016 


±0.019 


±0.002 


±0.003 


±0.000 


Rv 


±0.3 


T0.009 


±0.051 


±0.058 


±0.030 


±0.063 


T0.010 


Rvi 


±0.07 


±0.008 


±0.017 


±0.009 


±0.011 


±0.005 


±0.000 


Rvk 


±0.12 


±0.009 


±0.031 


±0.023 


±0.016 


±0.021 


T0.005 


Ar 


±0.008 


±0.002 


T0.011 


T0.015 


T0.004 


=F0.013 


±0.002 


A A. 


±0.007 


±0.005 


±0.017 


±0.013 


±0.011 


±0.012 


T0.002 


A(fi-V) 


±0.002 


T0.004 


±0.012 


±0.016 


±0.011 


±0.016 


±0.000 


A(V-7)c 


±0.006 


=F0.007 


T0.014 


T0.007 


T0.009 


=F0.007 


±0.002 


Total 




±0.080 


±0.079 


±0.074 


±0.051 


±0.072 


±0.012 


Fitting 




±0.023 


±0.058 


±0.053 


±0.038 


±0.020 


±0.007 



NOTE. — Estimated from JTH01 photometry. 
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rors in the reddening laws. Wegner (2003) determined Ry 
towards about 600 OB stars, from which we estimated an un- 
certainty of ±0.3 from the MAD of the distribution. We es- 
timated the dependence of R VI and R VK on the choice of R v 
from the extinction law of Cardelli et al. (1989), and found 
ARvi/ARv ~ 0.22 and AR VK /AR V ~ 0.84 for B-type stars. 
We adopted ±0.07 for the error in Ryi from the difference 
between our adopted value and Ryi = 124 from Dean et al. 
(1978). Similarly, we adopted ±0.12 as the error in R VK 
by considering our value with those advocated by Schultz & 
Wiemer (1975) and by Sneden et al. (1978). These are signif- 
icant contributors to the error budget for the highly-reddened 
cluster (NGC 25 16), but are much less significant for the other 
clusters. 

The next four rows list the errors resulting from zero point 
errors in the photometry. In this calculation, we assumed 
that the fundamental observed quantities are V, K s , B-V, 
and V—Ic (e.g., Stetson et al. 2003). The adopted values of 
the errors are shown in the second column. The error in K s 
was taken as the calibration uncertainty (0.007 mag) that was 
specified in the explanatory supplement to the 2MASS All 
Sky Data Release. 12 For the Pleiades and Praesepe, the errors 
come from intercomparisons between limited subsets of stars 
in common among various studies and from a consideration 
of the entries in Tables 2 and 3. For M67 and NGC 2516, we 
assumed the magnitude and color errors to be half the differ- 
ence between the two studies (Table 2). 

Besides the systematic errors in the photometry, distance 
estimation can be affected by the variability of pre-MS and 
MS stars (e.g., Scholz & Eisloffel 2004), which may be 
caused by rotational modulation of cool and hot regions on 
stellar surfaces. This issue can be addressed with multiple ob- 
servations for each star, but those kinds of data are not usually 
available. 

Errors in the metallicity and the helium abundance are im- 
portant contributors to the error in the distance (m-M)o,s- As 
shown in Table 5, individual spectroscopic studies generally 
have <~ 0.03 dex errors in the mean metal abundance, and this 
is propagated into the error of <~ 0.03 mag in distance modu- 
lus since A(m-M) /A[Fe/H] w 1. In M67, the error in the 
metallicity has no greater effect than the other systematic er- 
rors due to the many independent metallicity estimates. This 
contrasts well with the NGC 2516 case where only two stars 
were measured from high-resolution spectroscopy. The zero 
point errors in the photometry dominate the systematic errors 
in the photometric metallicity and the distance derived purely 
from the isochrones, independently of the spectroscopy. How- 
ever, our analysis demonstrates that the photometric metallic- 
ity and the distance can be derived with comparable accuracy 
to those from the spectroscopic studies. In all cases, the fitting 
errors are smaller than the total systematic errors by factors of 
2 - 3 as shown in the tables. 

The cluster richness can also affect the distance determi- 
nation. If there are N genuine single cluster members with 
random photometric errors in colors of e, the error in distance 
would be approximately se(3N)~ - 5 ~ 2eN~° 5 , where s is a 
typical slope of MS. This error was already included in the fit- 
ting error, but its contribution is small for clusters with good 
photometry and N ~ 100. However, information on single 
cluster members is only available for a handful of open clus- 
ters, and the effects of binaries and foreground/background 
stars should be taken into account in distance estimation, as is 

12 See http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html. 



discussed in the next two sections. 

6.2. Cluster Binaries 

The filtering algorithm (§ 4.1) removes cluster binaries if 
they are sufficiently far from the MS. However, low mass ra- 
tio binaries would remain and might systematically reduce the 
estimated distance since they are brighter and redder than sin- 
gle cluster members. 

To determine the size of the bias in distance estimation and 
to evaluate the effectiveness of the filtering algorithm, we per- 
formed artificial cluster CMD tests using a solar-metallicity 
isochrone for an age of 550 Myr. We constructed each set of 
CMDs by generating 200 single stars or unresolved binaries 
in the color range used for the MS fitting, 0.4 < (B-V)q < 1 .3. 
Primaries and secondaries in the binaries, as well as single MS 
stars, were randomly drawn from the mass function (MF) for 
M35 (Barrado y Navascues et al. 2001). The lower mass limit 
for the secondaries was set to be the hydrogen-burning limit at 
0.08M Q . Each simulation is characterized by the binary frac- 
tion (BF), which is defined as the number of binaries divided 
by the total number of systems in the above fitting range. For 
example, at BF = 50%, one-half of the data points correspond 
to single stars and the other half represent unresolved binaries. 
We neglected multiple systems other than binaries due to the 
observed low frequency among solar-type stars (Duquennoy 
& Mayor 1991). We assigned equal photometric errors of 0.01 
mag in V, B-V, V—Ic, and V -K s and displaced them from 
the isochrone assuming a normal distribution. The artificial 
CMDs were then filtered, and the distance was derived in the 
usual way. 

Figure 15 summarizes the results of these simulations. The 
top three panels show the bias in the distance modulus as a 
function of BF. The thick solid line shows the median of the 
bias from 200 artificial cluster CMDs computed at each BF 
with intervals of 5%. Thin lines on either side are the first and 
third quartiles of these distributions. The dashed line indicates 
A(m-M)o = 0. The vertical error bar displays the quartile 
ranges for a BF of 50%, but with errors for each star of 0.03 
mag. 

The bottom panel of Figure 1 5 displays the ratio of the num- 
ber of data points remaining after filtering compared to the in- 
put number of single stars in the simulation. The lines have 
the same meaning as in the other panels. For example, at BF 
= 50%, the simulation yields /V nt /N s - mg i e ~ 1.3 after filtering, 
indicating that ~ 70% of binaries were eliminated by the fil- 
tering algorithm. 

As seen in Figure 15, the distance from the (B-V,V) CMD 
is less affected by the remaining binaries than those from 
(V -Ic,V) and (V -K S ,V). Since stars were rejected if they 
were filtered at least in one of the CMDs, the total number of 
remaining binaries is the same in the three CMDs with a bi- 
nary detection limit set by a certain binary mass ratio 13 (<~ 0.4 
in our simulation). However, binaries stand out more promi- 
nently in (V-I c , V) or (V-K s , V) than in (B-V, V) because the 
addition of the cooler secondaries results in a smaller change 
in the combined color in B-V than in the other two color in- 
dices. As also seen in Figure 15, larger photometric errors 
result in a bigger bias in distance estimation since binaries 
are more difficult to identify when the photometric errors are 
larger. 

The effects of binaries on the determination of the photo- 

13 Binary mass ratio is defined as the mass of the secondary (less massive) 
star divided by the mass of the primary star. 
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FIG. 15. — Effects of remaining binaries after photometric filtering from 
artificial cluster CMD tests (see text for details). Top three panels: Bias in 
distance modulus as a function of binary fraction in each CMD with assumed 
photometric errors of 0.01 mag. Bottom: Number ratio of all data points used 
in the fit to the input single stars. The thick solid line shows the median of 
these values from 200 artificial cluster CMDs at each binary fraction with 
intervals of 5%, and thin lines on either side are the first and third quartiles. 
The dashed line indicates A(m-M)o = 0. The differences are in the sense of 
shorter distances and more remaining binaries at a higher binary fraction. The 
vertical error bar displays the quartile ranges for a binary fraction of 50%, but 
with errors for each star of 0.03 mag. 

metric metallicity and reddening are displayed in Figure 16. 
The top to bottom panels display the biases in [M/H]e, 
[M/H] % 2, and E(B-V) x i, respectively. Error bars at BF = 
50% show the results with the photometric errors of 0.03 mag. 
As the BF increases, the MS fitting gives lower metallicity 
estimates. This is a direct result of the systematic errors in 
distance as seen in Figure 15. Recall that if an isochrone at 
a particular metallicity yields a distance from (B - V, V) that 
is longer than from (V-/ c , V) or (V-K s , V), the photometric 
metallicity is smaller than that of the isochrone. Increasing the 
BF leaves the distances from (B—V, V) relatively unaltered but 
gives shorter distances in the other two colors, thus yielding 
the lower photometric metallicity. Nevertheless, the amount 
of a bias in the photometric metallicity is A[M/H] X 2 w -0.03 
at BF = 50%. In comparison, we conducted some experiments 
in which the filtering was not performed and found that this 
reduced the metallicity estimates by A[M/H] X 2 w -0.30. The 
results from the bottom panel show that the reddening deter- 
mination is almost unaffected by the presence of unfiltered 
binaries. 

In Figure 16 a small offset of A[M/H] £ w 0.02 at BF = 0% 
shows the systematic error that was already included in the 
"calibration" error in [M/H] £ (Tables 9 - 12). This happens 
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FIG. 16.— Same as Fig. 15, but for the bias in [M/H] E , [M/H] x 2, and 
E(B-V) X 2 (top to bottom). The differences are in the sense of a lower metal- 
licity and a higher reddening estimate at a higher binary fraction. 



because we assumed a linear relation between metallicity and 
distance over a wide range of metallicity (e.g., Fig. 5); higher 
order terms in the relation may be required to eliminate this 
systematic offset. 

Table 13 summarizes the predicted bias in the distance, pho- 
tometric metallicity, and reddening induced by unfiltered bi- 
naries at BF = 50%, which is typical of open clusters (Pa- 
tience et al. 2002). We compare results by adopting different 
MFs for secondaries, as listed in the first column. We further 
tested several combinations of the various MFs including the 
Kroupa et al. (1993) MF for both primaries and secondaries. 
However, we found that the detailed form of the MF has no 
greater impact on the distance determination than the inter- 
nal precision of photometry. The adopted random errors in 
the photometry are shown in the second column, while the 
remaining columns display the systematic errors. 

In Praesepe, many solar-type stars have been observed for 
binarity from IR speckle, direct imaging, and spectroscopic 
studies, covering a wide range of orbital periods (Mermilliod 
& Mayor 1999; Bouvier et al. 2001; Patience et al. 2002). 
As a result, its MS (Figs. 4 and 6) is fairly clean, suggesting 
that our binarity information in this cluster is almost com- 
plete at least for solar-type stars. In 0.4 < (B-V) < 0.8, 
there are 46 known binaries out of 1 16 systems, yielding a BF 
of 40%. At this value, our simulation results in Figures 15 
and 16 yield A(m-M) os = -0.006, A[M/H] £ = -0.014, 
A[M/H] X 2 = -0.037, and AE(B-V) = +0.003. These val- 
ues are consistent with those from the comparison of two sub- 
samples, computed before and after excluding known cluster 
binaries (Tables 7 and 8). 

The Pleiades has also been extensively surveyed on binarity 
(e.g., Mermilliod et al. 1992; Bouvier et al. 1997). However, 
there are many stars remaining above the MS even after ex- 
cluding all of the known binaries, suggesting that many bi- 
naries were probably left undetected by the previous surveys 
(see Fig. 7). In NGC 2516, there exists limited information on 
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TABLE 13 

Bias Induced by Unfiltered Binaries from Artificial Cluster 
TestatBF=50% 



MF a 


& phot 


A(m-M)o 


A[M/H] £ 


A[M/H] x2 


AE(B-V) X 2 


M35 


±0.01 


-0.010 


-0.007 


-0.033 


+0.002 


M35 


±0.03 


-0.028 


-0.034 


-0.050 


+0.003 


flat 


±0.01 


-0.006 


+0.004 


-0.020 


+0.002 


Salpeter b 


±0.01 


-0.010 


-0.006 


-0.023 


+0.000 



"Stellar mass function for secondaries. Primaries were generated from M35 
mass function (Barrado y Navascues et al. 2001). 
b Salpeter (1955). 



binarity only for B- and A-type stars from the spectroscopic 
studies (Abt & Levy 1972; Gonzalez & Lapasset 2000). 

To estimate the BF for each cluster, we first counted the to- 
tal number of stars that were 0.3 - 1.0 mag brighter in V than 
the cluster MS in (B-V,V). For NGC 2516, we subtracted 
field foreground/background stars based on the distribution of 
(m-M)o (see next section) from the JTH01 catalog before 
counting likely cluster binaries. Within the above AV range, 
the observed BFs are 19 ± 5% for the Pleiades and 18 ± 2% 
for NGC 25 16, while that of Praesepe is 19 ± 4%. If we com- 
bine these estimates with the total observed BF of 40% for 
Praesepe, the total BFs for the Pleiades and NGC 2516 would 
be 38 ± 1 1% and 38 ±4%, respectively. The errors represent 
those from counting statistics. 

Previous studies have also noted that Praesepe and the 
Pleiades have the similar BF among G- and K-type stars over 
a certain range of orbital periods (Mermilliod & Mayor 1999; 
Bouvieret al. 2001; Patience et al. 2002). However, the above 
BF estimates based on the observed fraction of stars with 
AV > 0.3 are uncertain for several reasons. For example, the 
distribution of binary mass ratio or the secondary MF can be 
different from one cluster to another. If we assume a flat MF 
for secondaries, we would derive BF s=s 50%. On the other 
hand, the M35 MF predicts quite a high BF (> 90%) and pro- 
duces too many binaries near MS (low mass ratio binaries) 
compared to the observed distribution of stars. It should also 
be noted that the BF in the simulation only concerns photo- 
metric binaries. In other words, blending of physically un- 
related stars could increase the BF, while binaries with large 
angular separations could reduce the BF. As a result, the BF 
depends on a number of cluster properties, as well as a spe- 
cific design of a photometric survey. In addition, a substan- 
tial fraction of binaries tend to have equal-mass companions 
(e.g., Halbwachs et al. 2003; Pinsonneault & Stanek 2006, and 
references therein). These binaries would be well-separated 
from the cluster MS, so they can be easily detected and re- 
moved by the filtering algorithm. The net effect would be an 
overestimation of the bias in distances and photometric metal- 
licities in the simulation for a given BF. 

From the above considerations, we adopted 30% - 50% as 
the range of ± 1 a formal errors in the BF for the Pleiades and 
NGC 2516. From the simulation result (0.01 mag error, M35 
MFs for both binary components), this BF yields A(m-M)s = 
-0.007 ±0.003, A[M/H] £ =-0.020 ±0.015, A[M/H] x2 = 
-0.036 ±0.013, and AE(B-V) X 2 = +0.003 ±0.001. These 
values should be subtracted from the MS-fitting results in Ta- 
bles 7 and 8 to correct the effects of binaries. Additional 
simulations showed that A[M/H] x2 = ±0.008 ± 0.004 and 
AE(B-V) x i =±0.003 ±0.001 in the case of the fixed clus- 
ter distance. Bias in distance at the photometric metallicity 



would be the quadrature sum of A(m-M) s and the error con- 
tributions from A[M/H] x2 (or A[M/H] £ ) and AE(B-V) x i. 

6.3. Field Star Contamination 

In clusters such as NGC 2516 (Fig. 9), a significant num- 
ber of foreground/background stars are present near the MS. 
Because the number density of these stars in each CMD in- 
creases towards fainter magnitudes, the distance from MS fit- 
ting can be systematically overestimated. 

In Figure 17 we show the distribution of individual dis- 
tance moduli for NGC 2516 from the JTH01 photometry in 
(B-V,V). Plotted are the number of stars per 0.1 mag in- 
terval in distance modulus. The same isochrone and color 
fitting ranges in Figure 9 were used. In the top panel, the his- 
togram represents the number of stars in the catalog before 
photometric filtering. It shows the density peak that repre- 
sents the cluster MS, the binary sequence that extends toward 
shorter distance from the MS, and the field star distribution 
that was fitted by an exponential function (solid line). The 
fit excluded stars near the MS and the cluster binaries, and 
was performed over a wider range in (m-M)o than shown 
here. In the middle panel the hatched histogram represents 
the distribution of the stars that remained after filtering from 
the catalog. This is compared with an open histogram, which 
was found by subtracting the exponential curve from the his- 
togram of the full data in the top panel. In the bottom panel we 
show the summed distribution of the RV members (Terndrup 
et al. 2002) and the Chandra X-ray-detected sources (Dami- 
ani et al. 2003) as an open histogram. The distribution of these 
stars after filtering is shown as a hatched histogram. 

One of the most conspicuous features in the middle panel 
of Figure 17 is the effective removal of cluster binaries from 
filtering. This can be seen from a deficit of data points in the 
hatched histogram against the open one at shorter distances 
from the MS. In addition, good matches are found between 
these two histograms near the MS, indicating that the filtering 
algorithm has worked correctly in the foreground/background 
subtraction. 

The RV and X-ray samples contain many cluster binaries 
(Fig. 17, bottom panel), but they are mostly free from field 
star contamination. On the other hand, the full catalog data 
(Fig. 17, top panel) are contaminated by both cluster bina- 
ries and foreground/background stars. Therefore, the effects 
of field star contamination after filtering can be estimated by 
comparing distances, metallicities, and reddening values de- 
rived from each filtered set of the data (Fig. 17, hatched his- 
tograms). From Tables 7 and 8, the weighted mean differ- 
ences in these quantities from the JTH01 and S02 photometry 
areA(m-M) s = ±0.010±0.002,A[M/H] £ =±0.020±0.003, 
A[M/H] x2 =±0.038 ±0.031, and AE(B-V) x i = -0.003 ± 
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FIG. 17. — Distribution of individual distance moduli for NGC 2516 from 
the JTH01 photometry in (B - V, V). Plotted are the number of stars per 
0.1 mag interval in distance modulus. The same isochrone and color fitting 
ranges in Fig. 9 were used. Top: Distribution of the full data in the cata- 
log. The solid line is an exponential fit to the histogram well beyond the MS 
and the binary sequence to model the field star distribution. Middle: The 
hatched histogram represents the distribution of the stars that remained after 
filtering. For comparison, the distribution of the full data set after subtract- 
ing backgrounds with the exponential curve in the top panel is shown as an 
open histogram. Bottom: The open histogram is a summed distribution of 
the radial velocity members (Terndrup et al. 2002) and the Chandra X-ray- 
detected sources (Damiani et al. 2003), and the distribution of these samples 
after filtering is shown as a hatched histogram. 

0.007. The errors are half of the differences in the estimated 
bias from the two photometric data sets. Within the estimated 
errors, these values are of the same size, but opposite direc- 
tion as the binary corrections at BF = 40%. In other words, 
these two biases cancel out each other. 

It is possible that the RV or X-ray samples are more heavily 
contaminated by cluster binaries than the full catalog data. 
However, we found that they have a similar fraction of stars 
(~ 19%) in 0.3 < AV < 1.0 mag above the MS. This implies 
that the RV and X-ray samples would have a similar total BF 
unless they have an excess of low mass ratio binaries. 

6.4. Comparison with Previous Estimates 

In Table 14 we summarize our metallicity and reddening 
estimates from the \ 2 minimization, and compare them with 
the literature values. The [M/H] and E(B-V) for the Pleiades 
are those derived at the average geometric distance in Table 1 . 
These estimates for the Pleiades were also corrected for the 
bias induced by the binary contamination. This correction was 
applied to the result from the data set that originally contained 
all known binaries. For NGC 2516, we did not apply the bias 
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FIG. 18. — Comparison of the photometric [M/H] x 2 and the spectroscopic 
[Fe/H] (top) and that of the photometrically determined E(B-V) X 2 and the 
literature value {bottom) from Table 14. The diagonal lines indicate identity. 

correction for both binaries and field star contamination for 
the reason explained in the previous section. The errors are 
the quadrature sum of the fitting and total systematic errors in 
Tables 9-12. For the Pleiades and NGC 2516, we addition- 
ally include the errors in the bias from binaries. Our estimates 
for M67 and NGC 2516 are the weighted averages from the 
two photometric data sets for each cluster. We adopted half 
of the difference in these estimates as the 1 a error in the 
case where they do not agree with each other within the esti- 
mated errors. In the third row of each quantity, weighted av- 
erages of our estimates and those in the literature are shown; 
however, the averages remain close to the literature values in 
most cases. Our estimated distances at these average metallic - 
ities and reddening values are shown and compared with the 
Hipparcos measurements (Robichon et al. 1999; van Leeuwen 
1999) in the last three rows. We discuss the Hipparcos mea- 
surements in § 7.2. 

In the top panel of Figure 18 we compare the photomet- 
ric metallicities with the spectroscopic values in Table 14. 
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TABLE 14 

Summary of Open Cluster Parameters 





Praesepe 


Pleiades 


M67 


NGC2516 


[Fe/H] 


This paper" 
Literature 
Weighted mean 


+0.20 ±0.04 
+0.11 ±0.03 
+0.14 ±0.02 


+0.04 ±0.05" 
+0.04 ±0.02 
+0.04 ±0.02 


-0.02±0.05 c 
+0.00 ±0.01 
+0.00 ±0.01 


-0.07 ±0.06" 
+0.01 ±0.07 
-0.04 ±0.05 


E(B-V) 


This paper" 
Literature 
Weighted mean 


-0.003 ±0.008 
0.007 ±0.002 
0.006 ±0.002 


0.022 ±0.014 b 
0.032 ±0.003 
0.032 ±0.003 


0.042 ±0.019 c 
0.041 ±0.004 
0.041 ±0.004 


0.125±0.014 d 
0.1 17 ±0.002 
0.1 17 ±0.002 


(m-M)o 


This paper 6 
Robichon et al. (1999) 
van Leeuwen (1999) 


6.33 ±0.04 
6.28±0.13 
6.37 ±0.15 


5.63 + 0.02 1 

5.36 ±0.06 

5.37 ±0.07 


9.61±0.03 c 


8.03±0.04 d 
7.70 ±0.16 



"Solution from the \ 2 minimization. 

'Solution at the average geometric distance (Table 1). 

c Average solution from MMJ93 and S04 photometry with single cluster membership by S04. 
d Average solution from JTH01 and S02 photometry. 
e Estimated at the weighted mean [Fe/H] and E(B-V). 
f Average geometric distance from Table 1. 



The comparison yields ([M/H] %2 - [Fe/H] spec ) = +0.01 ± 
0.03 (s.e.m.). Although our photometric metallicity for Prae- 
sepe is <~ 2 a higher than the spectroscopic value, the x 2 with 
respect to the identity line is 4.3 from these four cluster points. 
In other words, the errors in [M/H] x2 were correctly esti- 
mated. We investigate a photometric metallicity scale over 
a wider range of [Fe/H] in the next paper of this series (D. M. 
Terndrup et al. 2007, in preparation). 

The E(B-V) estimates from the x 2 minimization are plot- 
ted against the literature values in the bottom panel of Fig- 
ure 18. Here the mean difference between the two measures 
is (E(B-V) X 2-E(B-V) U ) =-0.006 ±0.012 (s.e.m.). The 
X 2 is 2.4, which again indicates that the errors were esti- 
mated correctly. However, our photometric reddening from 
BVIcK s data is not completely independent from the litera- 
ture values from UBV data since most of the literature esti- 
mates in Table 6 are from the similar methodology that uses 
the observed colors of the unreddened MS stars. We did not 
include U -band measurement in the isochrone calibration be- 
cause of its low reproducibility and relatively large systematic 
errors among different studies (Bessell 2005, and references 
therein). These properties make it difficult to reliably use the 
J/ -band data in the high-precision distance estimates as in this 
paper. 

7. OTHER ISSUES 
7.1. Are the K Dwarfs in NGC 2516 Blue, Too? 

In Figure 19 we display the MS of NGC 2516 in (B-V,V), 
which at an age of 140 Myr may exhibit the phenomenon of 
blue K dwarfs as seen in the Pleiades CMD (Stauffer et al. 
2003; see also Fig. 7). The left panel displays the photometry 
of JTH01 in (fi-V, V) while the right panel shows that of S02. 
In each panel, the open circles are those stars that remained 
after photometric filtering. The filled circles and plus signs 
show the stars that are RV members (Terndrup et al. 2002) 
or X-ray-detected sources (Damiani et al. 2003) if they were 
accepted or rejected, respectively, by the filtering. 

The solid lines in Figure 19 are 140 Myr isochrones with 
[Fe/H] = +0.01 and E(B-V) = 0.117, fitted over 0.40 < 
(B-V)o < 0.80. The isochrone matches the JTH01 photom- 
etry down to at least V « 17, but the SBBL02 photometry is 
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B - V B - V 

FIG. 19. — Two different interpretations on young K dwarfs in NGC 2516. 
Photometry of JTH01 (left) and that of S02 (right) are shown. Open circles 
are the stars that remained after photometric filtering from the full catalogs. 
Filled circles and plus signs show the stars that are RV members (Terndrup 
et al. 2002) or X-ray-detected sources (Damiani et al. 2003) if they were 
accepted or rejected, respectively, by the filtering on these samples. The same 
isochrone and color fitting ranges in Fig. 9 were used. Note that stars with 
the S02 photometry are bluer than the isochrone forB-V > 1.0. 

bluer than the isochrone forB- V > 1.0 (My > 6.5). This sys- 
tematic difference was previously noted when we compared 
the photometry from the two studies (Fig. 2). 

The individual distance modulus for each star is shown in 
Figure 20, plotted against (B-V)q, along with the Pleiades 
and Praesepe stars that have similar colors. The solid line in 
each panel is the median value from fits over 0.4 < (B-V)q < 
0.8 (note that our analysis above used the full color range 
for Praesepe). While the JTH01 data show a flat distribu- 
tion over 0.4 < (B-V)q < 1.3 as in Praesepe, the red stars 
in the S02 data systematically trend toward longer distances 
for (B-V)o > 0.8 where the Pleiades K dwarfs also depart 
from the isochrone. 
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FIG. 20. — Individual distance moduli of stars in the Pleiades (top), 
NGC 2516 (middle), and Praesepe (bottom) from (B— V, V) with respect to the 
isochrones with spectroscopic metallicities and reddening values in the liter- 
ature (Tables 5, 6). Only stars that remained after filtering are shown. Hori- 
zontal lines indicate median distance moduli derived in 0.4 < (B-V) < 0.8. 
While the JTH01 data show a flat distribution, the S02 data exhibit longer 
distances at redder colors where the Pleiades K dwarfs also depart from the 
isochrone. 



From these plots alone it is not possible to know whether 
the blue K dwarf phenomenon exists in NGC 2516 since the 
interpretation depends on which data set is used. Stauffer 
et al. (2003) argued that the departure from the ordinary MS 
in the Pleiades is connected with stellar rotation rate since 
the stars with the fastest rotation show the most anomalous 
colors (see their Fig. 15). In Figure 21 we display the dif- 
ference in (m-M)o for stars in NGC 2516 from the average 
distance modulus in each CMD against the projected rotation 
velocity. The filled circles are from the JTH01 photometry, 
while the open circles are from the S02. Positive A(m-M)o 
indicates that stars are fainter (or bluer) than the isochrone, 
and the dashed line shows A(m-M) = 0. The vsim data 
are only available up to V w 15 or (B-V)o < 1.0 (Terndrup 
et al. 2002), and only stars with (B — V)o > 0.8 are shown. 
Within this restricted color range, a few stars with rapid rota- 
tion in (B— V,V) are about 0.1 - 0.2 mag fainter than the slow 
rotators. By comparison, rapidly rotating stars with similar 
colors in the Pleiades are about 0.5 mag below the MS. This 
is suggestive that the stars in NGC 2516 may also exhibit the 
blue K dwarf phenomenon although the degree of departure is 
smaller than in the Pleiades. The opposite correlation is seen 
in (V-K S ,V) as also noted by Stauffer et al. (2003), and we 
see no significant correlation in (V-Ic, V). 
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FIG. 21. — Dependence of the photometric anomalies of K dwarfs in 
NGC 25 16 on the projected rotational velocity. The filled circles are from the 
JTH01 photometry, and the open circles are from S02 for 0.8 < (B — V)q < 
1.0. Displayed are the individual distance moduli with respect to the average 
distance in each CMD, and the dashed line shows A(m-M)o = 0. The same 
isochrone and color fitting ranges in Fig. 9 were used. Positive A(m-M)o 
indicates that stars are fainter than the isochrone. 

7.2. Hipparcos Distances 

In Figure 22 we compare our derived distance moduli with 
the Hipparcos parallax results from Robichon et al. (1999) for 
the Pleiades, Praesepe, and NGC 2516 (M67 was too distant 
for Hipparcos to measure its parallax). In the top panel the 
open circles show the comparison with (m-M)o.s at the spec- 
troscopic metallicity. The errors are the quadrature sum of the 
fitting and total systematic errors in Tables 9-12. For the 
Pleiades and NGC 2516, we additionally included the error in 
the bias correction for binaries. We also plot as an open trian- 
gle the comparison of (m-M)o.s and the geometric distance to 
the Pleiades from Table 1 and connect the two Pleiades points 
with a dashed line. The bottom panel shows the difference be- 
tween the two distance measurements. The filled circles rep- 
resent the comparison with (m-M) x 2 and are shown without 
error bars. 

While our distance to Praesepe shows a good agreement 
with the Hipparcos value, the 1999 Hipparcos distance to the 
Pleiades exhibits more than a ~ 3 a difference at either the 
spectroscopic or photometric metallicity even if we take into 
account all systematic errors. Similarly, NGC 2516 shows a 
<~ 2 a difference. These results extend the previous discus- 
sion by Pinsonneault et al. (1998) and Terndrup et al. (2002), 
but the current results are based on isochrones with improved 
color calibrations over wider color ranges and with additional 
(V-K s ) colors. 

In fact, Pinsonneault et al. (2000) derived a consistent dis- 
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FIG. 22. — Comparison of the Hipparcos parallax results (Robichon 
et al. 1999) and our derived distance moduli. Open circles are those for the 
Pleiades, Praesepe, and NGC 2516 (left to right) when distances are deter- 
mined at the spectroscopic metallicities. The open triangles are the compari- 
son with the geometric distance to the Pleiades (Table 1), and the two Pleiades 
points are connected with a dashed line. The errors are the quadrature sum 
of the fitting and total systematic errors in Tables 9-12 and those in the bias 
induced by binaries. Filled symbols are the comparison with (m—M) x z and 
are shown without error bars. The diagonal line indicates identity in the top 
panel. 

tance to NGC 2516 with the Hipparcos value, but with a con- 
siderably lower photometric metallicity ([M/H] = -0.26) than 
presented in this paper. This happened because B — V and 
V — Ic colors of their isochrones from the Yale color calibra- 
tion (Green 1988) were affected by large systematic errors in 
the color- r e ff relations. The photometric metallicity was re- 
determined to be [M/H] = -0.05 by Terndrup et al. (2002), 
who employed the same colors, but with a preliminary cali- 
bration of the isochrones that later appeared in our Paper II. 
As in this paper, they found a distance to NGC 2516 that was 
longer than the Hipparcos value. 

In the Pleiades, we achieved an excellent fit in all CMDs 
and derived a photometric metallicity and distance that are 
consistent with the spectroscopic determinations and the ge- 
ometric distances from a variety of methods. In our analy- 
sis, we excluded B — V colors where the cluster K dwarfs are 
anomalously blue. Percival et al. (2003) presented BV(RI)c 
photometry for a sample of local stars with accurate Hip- 
parcos parallaxes and developed a model-independent MS- 
fitting method based on the relation between luminosity and 
metallicity for these stars. In order to select unevolved stars, 
however, they picked targets with 0.7 < (B — V) < 1.0 and 
found a Pleiades distance in (B - V, V) that was longer than 
that from (V-Ic,V), indicating a photometric metallicity of 
[M/H] w -0.3. This led them to suggest that the spectro- 
scopic metallicity of the Pleiades was in error, even though 
it is well established (Table 5). As they subsequently realized 



(Percival et al. 2005), their analysis was seriously affected by 
the anomalous B—V colors of the K dwarfs. In the later study, 
they showed that distances derived with V—Ic and V-K s were 
consistent with the correct distance for the Pleiades. 

The disagreement with the Hipparcos values for the 
Pleiades and NGC 2516 has been attributed to a local zero 
point error of the individual stellar parallaxes in the Hippar- 
cos catalog (Pinsonneault et al. 1998; van Leeuwen 2005). A 
more elaborate reduction of the Hipparcos data (van Leeuwen 
& Fantino 2005) would help to resolve the difference from the 
MS-fitting distances. 

7.3. The Helium Abundance of the Pleiades 

Because the luminosity of the model isochrones at fixed T e s 
is sensitive to the helium abundance (dM\, \/dY w 3), we can 
derive the helium abundance if the cluster distance and metal- 
licity are known. The isochrones employed here assume that 
the Sun's metal abundance is Z Q = 0.0176 and that the solar 
helium abundance in models lacking diffusion is F Q = 0.266. 
This yields AT/AZ =1.2 with the primordial helium abun- 
dance of Y p = 0.245 ± 0.002 (see Paper I). 14 We chose the so- 
lar models with no diffusion because microscopic diffusion 
would not be important for young clusters such as the Hyades 
and the Pleiades. In Paper I we derived the Hyades helium 
abundance, Trades = 0.271 ±0.006, using the Hipparcos dis- 
tance to the Hyades (de Bruijne et al. 2001), combined with 
the metallicity estimate from Paulson et al. (2003) and the 
masses of the components of the eclipsing binary vB 22 (Tor- 
res & Ribas 2002). Since the standard solar model including 
the effects of helium diffusion and rotational mixing implies 
an initial solar helium abundance of T0,; n it = 0.274 (Bahcall 
et al. 2001), the Hyades abundance is equal within the errors 
to the Sun's value at the time of formation, even though the 
Hyades has ^35% higher metal content than the Sun. 

At the spectroscopic metallicity for the Pleiades, [Fe/H] = 
+0.04 ±0.03, we derived {m-M\ s = 5.664 ±0.041, where 
the error is the quadrature sum of the fitting and total sys- 
tematic errors from Table 10, excluding the error contribu- 
tion from helium. It also includes the error in the bias cor- 
rection for binaries. Our distance is longer than the geomet- 
ric distance, (m-M) = 5.632 ±0.017 (Table 1), which in- 
dicates that the helium abundance of the Pleiades is higher 
than the assumed value in the model. Since the isochrones 
have Y = 0.268 at the metallicity of the Pleiades, the differ- 
ence between the two distances yields Y = 0.279 ±0.015 for 
the Pleiades. 

The Pleiades helium abundance is the same as both the ini- 
tial solar and the Hyades' abundances within the errors. To 
know whether there is a dispersion in the local helium abun- 
dance, it would be necessary to reduce the size of the system- 
atic errors in the MS fitting either by improving the photo- 
metric calibration or by obtaining more accurate cluster metal 
abundances. However, we can certainly rule out extreme he- 
lium abundances (Y w 0.34) that had been advocated to ex- 
plain the shorter Hipparcos distance to the Pleiades (Belikov 
et al. 1998). From the primordial helium abundance, the ini- 
tial solar value, and the helium abundances for the Hyades and 
the Pleiades, we derive AY / AZ = 1.42 ±0.17. 

When future astrometric missions such as Gaia (Perryman 
et al. 2001) provide highly accurate distances, the error in 
the helium abundance will be predominantly determined by 

14 We assumed that the helium abundance is a function of the heavy- 
element content Z as given by Y = Y P + (AY/AZ)Z. 
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the accuracy of the MS fitting. In such cases, cosmic scat- 
ter of AF = 0.02 could be detectable with the current ac- 
curacy of the MS-fitting technique. Furthermore, the error 
in the helium enrichment parameter scales as e(AF/AZ) w 
0.3(A^/100)-°' 5 (cr AJ /0.05 mag), where N is the total number of 
clusters and er M is the error in the distance from the MS fit- 
ting. A more accurate value of AF/AZ requires that many 
input systematic errors in the MS fitting be kept under con- 
trol. Stellar abundances in particular are susceptible to sys- 
tematic errors from the adopted r e ff scales. This issue will be 
addressed in the next paper of this series (D. M. Terndrup et 
al. 2007, in preparation). 

8. SUMMARY AND DISCUSSION 

We have generated a set of isochrones near solar metallic- 
ity (-0.3 < [Fe/H] < +0.2 in 0.1 dex increments) and have 
calibrated them by correcting the Lejeune et al. (1997, 1998) 
color- r e ff relations to give a precise match to the Hyades MS 
at [Fe/H] = +0.13 ±0.01 and (m-M) Q = 3.33 ± 0.01 (Pa- 
per II). We tested the hypothesis that the empirical color- T e g 
corrections improve the accuracy of distance estimation using 
the four well-studied open clusters - Praesepe, the Pleiades, 
M67, and NGC 2516. 

Either by adopting cluster metallicities and reddening val- 
ues in the literature or by constraining them from the multi- 
color CMDs, we showed that the calibrated set of isochrones 
not only matches the cluster MS precisely but also gives in- 
ternally consistent distances from (B-V,V), (V-Ic,V), and 
(V-K S ,V) CMDs. The internal precision of our derived dis- 
tances is about a factor of 5 better than the case without the 
color calibrations. The current result supports our previous 
assertion (Paper II) that the largest source of isochrone mis- 
matches stems mainly from the conversion of theoretically 
predicted quantities into magnitudes and colors. 

We also found that our photometric metallicities are in good 
agreement with spectroscopic estimates although the largest 
difference was found for Praesepe, where the photometric 
metallicity was <~ 2 a higher than the spectroscopic value. In 
addition, we found an excellent agreement for the Pleiades be- 
tween our derived distance and the geometric determinations 
other than the original Hipparcos measurements, confirming 
the validity of the empirically calibrated isochrones. 

We explored various systematic errors in distance, photo- 
metric metallicity, and reddening estimates. The distance 
from MS fitting is sensitive to the metallicity and helium 
abundance but can be estimated to an accuracy of ± 1 % (fitting 
error) and ±2% (systematic error) if there exists a good spec- 
troscopic abundance measurement. When we determine clus- 
ter metallicity and/or reddening with our isochrones alone, the 
distance becomes sensitive to the size of the systematic zero- 
point errors in the photometry. With reasonable sizes of the 
zero-point errors, however, distances can be determined with 
comparable accuracy to those at the spectroscopic metallici- 
ties. 

We showed that the Pleiades K dwarfs are bluer than the 
calibrated isochrone in (B-V,V), confirming the previous 
study by Stauffer et al. (2003). We further showed that the 



red stars in NGC 2516 are bluer than the isochrone when the 
S02 data are used; however, the JTH01 data closely match 
the isochrone down to the limit of the isochrone calibration. 
Using a limited set of vsin/ data, we found that a few stars 
with high vsin/ are fainter by 0.1 - 0.2 mag than the slow 
rotators in 0.8 < (fi-V)o < 1.0. This suggests that the K 
dwarfs in NGC 2516 may also exhibit the blue K dwarf phe- 
nomenon, although the degree of departure is smaller than in 
the Pleiades. More vsin / data for lower mass stars and careful 
photometric standardization are required to set the upper age 
limit of the anomalous colors of active K dwarfs. For young 
open clusters, careful attention should be given to this kind of 
anomaly to avoid erroneous results in distance. 

We also derived a helium abundance for the Pleiades of Y = 
0.279 ± 0.015, based on the average geometric distance to the 
Pleiades (Table 1). This value is equal within the errors to the 
Sun's initial helium abundance and that of the Hyades. This 
result implies that there is no observed cosmic scatter in the 
helium abundance at the level of the accuracy in the current 
estimation. To decrease the size of the error, more accurate 
abundances for the cluster and smaller systematic errors in 
the photometry are required. 

There are distinct limits to the current calibration. The de- 
pendence of the MS luminosity and colors on metallicity can- 
not be inferred from one calibrating cluster alone and must 
be either assumed to be correct in the models (as in this pa- 
per) or else checked against local field stars with known dis- 
tances and metallicity (D. M. Terndrup et al. 2007, in prepara- 
tion). The Hipparcos data in the Hyades are only complete to 
V w 1 1 (M v w 8), so the empirical calibration in Paper II ex- 
tends only down to ~ O.7M . In addition, the Hyades has an 
age of about 550 Myr, so younger clusters should be used to 
calibrate the upper MS for stars more massive than the Hyades 
turnoff stars. All of these issues will be addressed with pho- 
tometry of low-mass members in Praesepe and a Persei in the 
later papers of this series. 
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APPENDIX 

THE HYADES PHOTOMETRIC CALIBRATION REVISITED 

In § 6.1 we estimated the systematic errors in distance, photometric metallicity, and reddening that would result from zero-point 
errors in the photometry of target clusters. Here we elaborate on the systematic errors in the Hyades photometry, which determine 
the accuracy of the calibration itself. 
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TABLE A 1 
Comparison of Johnson and 
Kron Transformations 



Cluster 


A(V-/) c a 


N b 


Hyades 


-0.007 ±0.007 


5 


Praesepe 


-0.009 ± 0.003 


7 


Pleiades 


-0.010 ±0.004 
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"Transformed Johnson minus 
Kron photometry from previous 
transformations employed in 
Paper II (see text). 

b Number of stars used in the com- 
parison. 
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FIG. A 1 . — Differences between the magnitudes and colors of the Hyades members used to calibrate the Hyades empirical isochrone (Paper II) and photometry 
of Joner et al. (2006). The differences are in the sense of our compiled data minus Joner et al. values. The solid line in each panel indicates identity. The Johnson 
V-I colors used in Paper II were transformed to the Cousins system using an updated transformation (see text). 

In our assembly of the Hyades photometry in Paper II, we employed the Bessell (1979) transformation between the Johnson 
and Cousins V-I systems, namely, 

(V-/)c = 0.778(V-/). 7 . (Al) 

This calibration has a small zero-point error, in that this produces colors that are somewhat too blue. Cousins (1980), for example, 
derived 

(V-7) c = 0.77(V-7)j + 0.01. (A2) 

As part of a detailed study of local field stars, we have assembled a large sample of stars with both Cousins and Johnson 
photometry (D. M. Terndrup et al. 2007, in preparation). These indicate 

(V-/) c = 0.779(V-/)j + 0.014. (A3) 

Although the size of the zero point is small, it produces a systematic error in distance of <~ 0.1 mag. This is because the Johnson 
data for the Hyades mostly have My < 6 where the average slope of the MS is AV/A(V-/) C w 7. 

As summarized in Table Al, there are a few stars in the Hyades, Praesepe, and the Pleiades that have V -I photometry on both 
the Johnson and Kron systems. These provide additional evidence that a small redward shift in the transformation of the Johnson 
colors is needed. The second column shows the mean magnitude difference in the sense of Johnson minus Kron after the former 
were placed on the Cousins system using the (incorrect) transformation in equation (Al). The Kron colors were transformed 
to the Cousins system using the cubic polynomial derived by Bessell & Weis (1987). Errors are the standard error of the mean 
difference, and the number of stars used in the comparison is shown in the last column. In all cases, the Johnson-to-Cousins 
transformation from equation (Al) employed in Paper II produces colors that are slightly too blue compared to those from the 
Kron system. 

The adoption of the new Johnson-to-Cousins transformation (eq. [A3]) modifies our Hyades calibration, making the Hyades 
empirical isochrone from Paper II slightly redder in the upper MS since the photometry for these stars is mainly on the Johnson 
system. The net effect is to make our color- T eS correction table (Table 2 in Paper II) redder by A(V-I) C = 0.012 for My < 6. 
The isochrones used in this paper included this correction, no adjustment for My > 8, and a linear transition between these two 
magnitude limits. We also employed equation (A3) in the assembly of the Praesepe and the Pleiades photometry in this paper. 

In Paper II, we did not include V—Ic from Taylor & Joner (1985) because we were unable to achieve a good transformation 
between their photometry and those on the Johnson system from other sources (Johnson et al. 1966; Mendoza 1967; Johnson et al. 
1968). The exclusion of their data prompted Taylor & Joner (2005) to discuss their photometry at length and to examine whether 
there were systematic errors in the source photometry used in our previous work. They first verified that their V—Ic values are 
about 0.02 mag redder than ours. Then, they identified two sources of the difference: the adjustment in the Johnson-Cousins 
transformation noted above (but derived independently here), and a systematic error in Mendoza's R—Ij photometry. 
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Recently, Joner et al. (2006) published new Cousins photometry of a large sample of Hyades stars. Their data set has the 
advantage over the previous collection in that it was taken on a single telescope and filter set. In Figure Al we plot the differences 
in the photometry for the stars used to calibrate the Hyades isochrone (Paper II) and the data in Joner et al. (2006). Compared 
to our values in Paper II, their photometry is brighter in V and slightly redder in B - V . In V - I c , the bright stars in our sample 
are bluer than the Joner et al. values, but both data are in agreement for the fainter stars. The mean differences and standard 
errors for the range used in the MS fitting of this paper (M v > 3.1) are AV =+0.026 ±0.002, A(fi-V) = -0.009 ± 0.001, and 
A(V-I) C = +0.004 ±0.002. 

These differences can be compared to those found previously (Table 2) between the Tycho photometry and our values in 
Paper II. Compared to the Tycho-1, our Hyades values are also fainter by 0.012 mag and are bluer in B-V by 0.009 mag. If 
we had adjusted zero points in V, B-V, and V-/ c of isochrones to the Joner et al. (2006) scales, the distances at spectroscopic 
metallicities would become longer than our estimated values in Table 7 by 0.01 mag for Praesepe, 0.03 mag for the Pleiades, and 
0.02 mag for both M67 and NGC 25 16. 
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